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Abstract. We renormalize six dimensional theory in the modified minimal subtraction (MS) 
scheme at four loops. From the resulting /^-function, anomalous dimension and mass anomalous 
dimension we compute four loop critical exponents relevant to the Lee-Yang edge singularity 
and percolation problems. Using resummation methods and information on the exponents of the 
relevant two dimensional conformal field theory we obtain estimates for exponents in dimensions 
3, 4 and 5 which are in reasonable agreement with other techniques for these two problems. The 
renormalization group functions for the more general theory with an 0{N) symmetry are also 
computed in order to obtain estimates of exponents at various fixed points in five dimensions. 
Included in this 0{N) analysis is the full evaluation of the mass operator mixing matrix of 
anomalous dimensions at four loops. We show that its eigen-exponents are in agreement with 
the mass exponents computed at 0{1/N‘^) in the non-perturbative large N expansion. 
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1 Introduction. 


Recently there has been renewed interest in analysing the scalar quantum field theory which is 
perturbatively renormalizable in six spacetime dimensions. This interest in primarily due in the 
main to the modern development of the original conformal bootstrap method, [DEIEIIIEIEIE], 
to study the fixed point structure of field theories in a nonperturbative way, (HusiiniE]. One 
of the aims is, for example, to ascertain whether a conformal window exists and if so for what 
range of parameters of the symmetry group. Such studies are not limited to two dimensions 
where the structure of the conformal group being infinite dimensional, is fundamentally different 
to the finite dimensional conformal group in d > 2 where d is the spacetime dimension. Instead 
models exhibiting conformal symmetry in three, four and higher dimensions are of interest. One 
motivation for such studies rested in part on applications to dualities in higher spin AdS/CFT’s 
[mils] as well as model building beyond the Standard Model. Indeed in terms of gauge theories 
the original study in quantum chromodynamics (QCD) of [H] suggested that for a certain range 
of the number of quarks there could be a nontrivial fixed point in strictly four dimensions. Such 
a fixed point, which is known as the Banks-Zaks fixed point, may have a connection with the 
chiral symmetry phase transition, M- Another motivation for studying conformal properties in 
higher dimensional theories rests in trying to generalize properties of theorems such as the two 
dimensional c-theorem, m, to analogues in three and four dimensions. Though to have some 
insight into such extensions one has to be aware of the fixed point structure of the underlying 
quantum field theory. In this respect there has been interest in tackling this problem in scalar 
(p^ theory in six and lower dimensions. See, for example, [16t flTl [T8| [19] for recent in depth 
studies of the 0{N) symmetric theory. For instance, the fixed point structure has been 
comprehensively studied perturbatively to three loops in [181120] . This has also subsequently 
been extended to the theory with an Sp{N) symmetry in |21| . 

Those works exploit the renormalization of the theory from a generation ago, [221 , 

when modern multiloop computational techniques were not available. One highlight of m 
[20] was the estimate of the conformal window. It is possible to compute order by order in 
perturbation theory the value of as a function of the spacetime dimension for which the 
stable infrared fixed point ceases to exist. This critical value is denoted by N^r- Evaluating the 
expression for Ncr in five dimensions the one loop result of 123 [IE] was that the leading order 
value is Ncr = 1038. Indeed the first examination of the conformal window in six dimensional 
0{N) (p^ theory was given in [25]. In [20] the three loop analysis reduced this to Ncr = 64. 
Such a large reduction gives credence to the hope that the value can be reduced further. Indeed 
there appears to be support for this in conformal bootstrap analyses, [261 [271 [28], and moreover 
the series for Ncr hr m appears to be convergent. However, partly in order to resolve this but 
mainly for applications to other problems, the primary aim of this article is to extend the work 
of [22[|23l[24] to four loops. This is now possible given the advances in methods to evaluate high 
loop massless Feynman diagrams and in particular 2-point functions in the period since [23[ 124] . 
Moreover, while this will involve a large amount of integration such a renormalization could 
not proceed in the absence of powerful symbolic manipulation languages as well as powerful 
computing resources. Therefore, we will not only provide a comprehensive analysis of the four 
loop structure of the renormalization group functions of cp^ theory in the modified minimal 
subtraction (MS) scheme, but we will also outline the computational algorithm used. In this 
respect we will provide the basic four loop massless 2-point function Feynman integrals which 
are central to the calculations to allow others to extend the programme to other models. 

In indicating our analysis will be extensive this not only means that we will focus on the 
0{N) symmetric theory but also simpler models where there is no or different symmetry. This is 
because of the connection of these theories to condensed matter or statistical physics problems. 
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The simplest such case is the theory with one scalar field. It is of interest because of the 
relation to the Lee-Yang edge singularity problem, [29]. When the coupling constant of the 
single field theory is purely imaginary, |29|, then the critical exponent o", which is determined 
at the d-dimensional Wilson-Fisher fixed point, is the main quantity of physical interest. It can 
be determined via the renormalization group functions of the theory. The exponent cr has been 
estimated by different techniques such as high temperature series and Monte Carlo methods for 
the discrete dimensions less than six. Several references to such work are [aolEIlESI, for instance. 
More recently the conformal bootstrap programme have been applied to it, |33|- With these 
evidently more powerful methods giving reasonable agreement it is therefore timely to extend 
the e-expansion analysis of critical exponents to O(e^) where we will use d = 6 — 2e throughout. 
Not only will this improve the estimates given in [231124| but we will follow the method of |32j 
where constrained Fade approximants were used. Central to this idea is the exploitation of the 
properties of the underlying two dimensional conformal field theory. There the corresponding 
critical exponents are known exactly. This information is used like a boundary condition on the 
Fade approximant and, as will be evident, will significantly improve exponent estimates for large 
values of e. It will transpire that this analytic perturbative approach will give exponent estimates 
which are not unreasonable in comparison with numerically intense alternatives. While these 
comments have been driven by the Lee-Yang edge singularity connection to theory, they will 
equally apply to a similar, in terms of the underlying quantum field theory, but different physical 
problem. This is the percolation problem which can be formulated in continuum </>^ theory by a 
multiplet of scalar fields with a group valued coupling constant, [3l|- In the appropriate replica 
limit the critical exponents at the Wilson-Fisher fixed point equate to the 1-state Fotts model 
and hence percolation. We will provide estimates for a large set of critical exponents in the 
dimensions between two and six. Again it will be the case that the estimates are in keeping with 
other approaches. 

That we draw attention to these applications in the same context as the recent 0{N) (j)^ 
theory studies is not unconnected. As outlined in [20| the 0{N) theory has to have a connectivity 
with an underlying two dimensional conformal field theory in respect of the various critical 
points which emerge. If ultimately the d-dimensional theory can be identified with such a two 
dimensional theory then it will open the possibility of extended constrained Fade approximants 
to the exponents of the various fixed points of the 0{N) theory. Though it is not currently clear 
what the application to a physical problem is at the moment. One other motivation in [18] was 
to make the connection of theory with 0{N) theory in the dimension range 4 < d < 6 at 
the Wilson-Fisher fixed point via the large N expansion in d-dimensions. It is well known that 
in 2 < d < 4 the 0{N) theory is in the same universality class as the 0{N) nonlinear a model 
and the three dimensional Heisenberg ferromagnet. Above four dimensions it turns out that at 
the Wilson-Fisher fixed point the d-dimensional theory is in the same universality class as one 
of the 0{N) fixed points. So using the large N expansion we can evaluate the four loop 
d-dimensional critical exponents and compare them with those known from the explicit large 
N expansion, [SaiMlEZllM]. This will provide a nontrivial check on our perturbative results. 
Indeed we will compute the full mixing matrix for the mass operators in the 0{N) case and show 
the subtlety in connecting the mass eigen-critical exponents with the two separate mass critical 
epxonents in the d-dimensional 0{N) p^ theory. In this context it is interesting to note that 
the determination of the exponent rj at [37] . was originally evaluated as a function of 

d by using a d-dimensional conformal bootstrap method but in an analytic as opposed to the 
modern numerical approach. 

The article is organized as follows. We summarize the relevant background to performing the 
four loop renormalization of the basic six dimensional p^ theory in section 2. The notation used 
throughout is given there as well as the technical details of how the computation was organized. 
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In particular we emphasise that all the basic renormalization constants can be extracted purely 
from an evaluation of the 2-point function. Results of this method are summarized in the 
subsequent section where the renormalization group functions are recorded for various symmetry 
configurations. Included in this are theories with SU{Nc) symmetry which were examined as 
early toy examples of the strong interactions. Sections 4 and 5 provide details respectively of 
the results for the Lee-Yang edge singularity and percolation problems. The final main problem 
we analyse is the full four loop structure of the 0{N) version of </>^ theory which is provided 
in section 6. Conclusions are given in section 7. In addition two appendices are provided. The 
first gives the e expansions of the relevant basic but nontrivial integrals we needed at four loops. 
The second extends the three loop results for various values of N given in the appendix A of 
[28] to four loops but also includes the mass eigen-critical exponents for each fixed point as an 
expansion in e. 


2 Background. 


In this section we outline the technical aspects of the four loop renormalization of six dimensional 
theory. We will consider the theory in different guises depending on the particular application 
required to extract, for instance, critical exponents. These will depend essentially on different 
decorations of the underlying Lagrangian, which is 

L = 1 (2.1) 

with various symmetry groups. Here g is the coupling constant which is dimensionless in six 
dimensions. For example, both the Lee-Yang edge singularity and percolation problems can be 
accommodated with the more general theory, [23] l24] , 

L = \ ( 2 . 2 ) 

where there are a multiplet of fields and a group theory tensor which is totally symmetric in 

its indices in addition to a coupling constant g. The Lagrangians for theories with more than 
one set of fields will be discussed in later sections but the underlying calculational procedure is 
effectively the same as that introduced here. For (12.21) we use the same assumptions as [23 E4] 
for the renormalizability of the Lagrangian. Briefly this reduces to two observations. First, for 
2-point self-energy graphs, including subgraphs, the product of two coupling tensors satisfies 

^ilhi4^i2i3U ^ ^ 2 ^ 1*2 ( 2 . 3 ) 


where we use a different and more systematic notation here compared to [231 [24] for what 
corresponds to the group theory Casimirs. The second assumption in [231124] is that the product 
of coupling tensors is such that if the product has three free indices then it is proportional to 
itself. In this way (12.21) will clearly be renormalizable. Given this it might be thought that there 
is a sizeable number of group invariants which can appear at high loop order. This is not the 
case as it transpires that to four loops the following Casimirs suffice to write the renormalization 
group functions in a concise form. These are 


fjihh ^ki2i3 _ ^* 1 * 3*5 ^* 2 * 4*6 _ 
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^**1*2 fjkillh2 ^*1*5*6 ^*2*7*8 ^*3*5*9 ^*447*10 ^*6*8*11 ^^*9*10*12 _ 

^**1*2 ^i*3*4^fc*6»12 ^*1*5*6 ^*2*7*8 (^*3*5*9 ^*4*7*10 ^*8*11*12 ^*9*10*11 _ 

^*n*2^i*3*4^fc*5*12 ^*1*5*6 ^*2*7*8 ^*3*9*12 ^*4*10*11 ^*6*7*10 ^*8*9*11 _ 

^*n*2^i*3*4^fe*8*12 ^*1*5*6 ^*2*7*8 ^*3*5*9 ^*4*10*11 ^*6*7*10^*9*11*12 _ 

^**1*2 ^^*3*4 ^fc*llil2 ^*1*3*5 ^*216*7^*4*6*8 ^*5*9*10 ^*7*9*11 ^*8*10*12 _ 

^** 1 * 2 ^^* 314 ^^ 15 * 12 ^^* 1 * 5 * 6 ^* 217 * 8 (^ 13 * 9112 ^ 14 * 7110 ^ 16110111 (^ 18 * 9*11 _ 

(^** 1 * 2 (^^* 3 * 4 (^^ 16 * 12 (^* 1 * 5 * 6 (^* 21718 (^ 13 * 519 (^* 417 * 10 (^ 1819*11 (^* 10 * 11*12 _ 

(^1*1*2(^^1314(^^*11112(^11*5*6(^*217*8(113*519(^1417*10(^16110111(^18*9*12 _ 


(2.4) 


where we have included the previous invariants up to three loops of [231 El]. So at most there 
are nine new group invariants at four loops in principle. Though it will turn out that while 
these nine correspond to different graph topologies the four loop Casimirs, Tg^, do not all have 
distinct values. The syntax here is that the first label on Tn corresponds to the number of 
tensors in the product and the second label distinguishes between different invariants. If 
the first label is even then it arises in a 2 -point function while an odd label indicates an invariant 
which will appear in the coupling constant renormalization. Labelling the actual group indices in 
each product in (I2.4|l with subscripts allows one to straightforwardly construct the underlying 
topologies graphically. Each paired index corresponds to a propagator and these will join at 
vertices defined by each coupling tensor. To make contact with the notation used in [231 El] for 
the invariants up to three loops, we note that T 2 = a, T 3 = /3, T 5 = 7 , T 71 = 6 and T 72 = A 
when one observes that the defining graphs given in Figure 1 of [23] appear in reversed mirror 
image. We prefer the notation T„ here to avoid confusion with the standard notation of various 
critical exponents which also use Greek letters. 

Although the renormalization of both (12.ip and (|2.2I) are the same we will focus the technical 
discussion on the former as the extension to tensor couplings is not onerous. For (12.1|) we have to 
determine the wave function and coupling constant renormalization constants. Once these have 
been established then the determination of the mass anomalous dimension follows as a corollary 
even though the Lagrangian is massless. Indeed we work with massless fields throughout as this 
renders the four loop /3-function accessible with a relatively minimal amount of computation. 
The first stage is to compute the 2-point and 3-point functions of (|2.1I) . The Feynman graphs are 
generated using the Qgraf package, [39]. There are respectively 1, 2, 10 and 64 one, two, three 
and four loops diagrams to evaluate to the simple pole in e in dimensional regularization which 
we use throughout. Indeed we note that the index labelling in (|2.4I) is based on the Qgraf 
output. We note that throughout our convention is d = 6 — 2e. For the 3-point vertex the 
respective numbers are 1, 7, 56 and 540 which is an order of magnitude increase at four loops 
compared to the 2-point case. However, we will compute the coupling constant renormalization 
constant for ( 12 .ip purely from the 2 -point graphs by exploiting certain properties of the specific 
field theory. It is based on the following observation. If for the moment one considers the massive 
extension of () 2 .ip then the propagator can be formally expanded in powers of m? where m is 
the mass via 

The first term on the right hand side of course corresponds to the massless theory. However, 
the second term represents the zero momentum insertion of the unit operator on a propagator. 
Diagrammatically for the self energy renormalization this corresponds to a 2-point function 
with a zero momentum insertion but more importantly this term would correspond to a 3- 
point function graph where one of the external legs has a nullified momentum. In other words 
this would be equivalent to a graph contributing to the coupling constant renormalization. As 
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we are computing in dimensional regularization and only interested in the MS renormalization 
scheme the coupling constant renormalization constant can be correctly extracted from this 
nullified external momentum configuration. Indeed in four dimensional gauge theories this is 
the standard procedure for three loop renormalization [iQnn]. One concern is that nullifying 
an external momentum in (j 2 . 1 jl could introduce unwanted infrared divergences which would be 
indistinguishable from ultraviolet ones in dimensional regularization. Indeed in four dimensions if 
a propagator was present in a massless graph which had the form of the second term of ( 12 . then 
it would be infrared singular. That it is not an issue for ()2.ip is because in six dimensions such 
a propagator is by contrast infrared safe even in scalar field theories unlike in four dimensions. 
Therefore, this observation radically reduces the number of Feynman graphs to be computed. 
One downside is that there are more integrals to determine but this is surmountable as we 
will indicate later. One concern which may arise with the expansion approach is that symmetry 
factors in the 2 -point function may not be set in such a way that after expansion the term relating 
to the vertex has the incorrect factor. We have checked that it is the case to three loops. We 
illustrate this at one loop. There the basic one loop self-energy graph has a symmetry factor of 
i but the one loop 3-point vertex has a unit symmetry factor. However, using (12.51) there are 
two O(m^) terms which cancel the factor of ^ so that everything tallies. One final point resides 
in the determination of the renormalization of the mass operator in (12.ip . It transpires that 
in this theory the mass and coupling constant renormalization are equivalent. This is implicit 
in the algorithm we discussed for the determination of the latter and has already been noted in, 
for instance, [23l[2l]. So in (12.ip no extra computation has to be carried out. However, for (12.2|) 
the mass and coupling constant renormalization constants are different. This is evident in the 
algorithm we have introduced. Both can be determined using our method of only computing 
the graphs for the 2-point function by extending (|2.5I1 and mapping each massless propagator 
with 

^ + "(pyr + (A:2)2 (2.6) 

after the graphs have been generated with Qgraf. We have included the group theory structure 
on the propagators and m? are parameters of the dimension of a mass. These parameters are 
included as a label to distinguish which term is which when one comes to extract the various 
wave function, mass and coupling constant renormalization constants in the sum of all the 
contributions to our 2-point function. That involving mi corresponds to the mass operator. 
The final term of ()2.6p with label m 2 corresponds to the zero momentum insertion for the 
renormalization of the coupling constant g. The index kg, in (12. 6 p is the index associated with 
the third leg of the 3-point function and is a fixed free external index. 

Having discussed the method to isolate the necessary graphs and integrals contributing 
to the wave function, coupling constant and mass operator renormalization constants we now 
discuss the method used to determine their divergent part. The approach is to use the Laporta 
algorithm, [52]. This is a systematic application of integration by parts to establish towers 
of relations between integrals which are then solved algebraically to express all integrals in 
terms of a relatively small set of basic integrals which are called masters. The values for these 
have to be determined directly as integration by parts can no longer be used. Once these are 
found then the integration algorithm is complete and we run an automatic determination of the 
respective Green’s functions to find the underlying renormalization constants. Crucial to this is 
the symbolic manipulation language Form, [43], and its threaded version Tform, [44], which we 
use extensively to handle the large amounts of algebra. This arises since the reduction to masters 
produces relations where the coefficients are large rational polynomials in d. Not only have these 
to be expanded in powers of e and the masters substituted but the denominators may contain 
factors of (d—6) which are termed spurious poles. Thus not only have the polynomial coefficients 
to be expanded to higher order in e but some masters need to be evaluated to a precision beyond 
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the simple pole. Even for the simple scalar theory, m, where there are no tensor integrals 
certain graphs become tedious to evaluate. For the present four loop computation we have used 
the Reduze package, [l5]. One useful feature is that the output relations between integrals 
can be readily converted to Form input notation and thus included as a module within the 
automatic Form evaluation. 

The final ingredient in the algorithm is the explicit values for the master integrals. This 
requires a method other than integration by parts and we have adapted several approaches 
which have been used for similar problems in other contexts such as |451146j . In addition we 
were able to exploit a feature of Reduze which means that we could do this systematically. 
By this we mean the following. When the relation between integrals are solved in Reduze 
the algorithm has some internal criterion for determining which integrals are to be used as the 
masters in terms of which all other integrals are expressed. This is not always the best choice 
for the calculation of interest. For instance, some of the masters may be simple self-energy 
graphs where each subgraph is itself a lower loop self-energy graph or products of lower loop 
graphs. Such integrals are straightforward to determine and are retained in the set of masters. 
However, more difficult integrals remain such as those which are primitively divergent and it is 
their simple and only pole which will be needed. The feature of Reduze which one exploits is 
that one can specify a set of integrals which the package identifies as the masters. This is done 
after the initial reduction has been determined and a database constructed, [46]. Therefore, we 
have chosen a set of basic master integrals which, aside from those which are simple to evaluate, 
are hnite using Weinberg’s theorem, m- This technique has been elaborated on recently in 
|48j . While this may seem to resolve the computational algorithm one has to be careful. This 
is because of the problem of spurious poles in e which means that not only the leading term of 
these finite integrals are required, but sometimes also several terms in powers of e. However, 
one also has to have some information on a master choice to be able to solve for these and other 
required coefficients. 

To do this we use the known four loop self-energy master integrals given in [36] . These were 
determined by application of the glue-and-cut method of five loop primitive massless vacuum 
diagrams to varying orders in the e expansion. The results are consistent with a subsequent 
independent numerical sector decomposition evaluation given in |49] . Analytic evaluations were 
also developed thereafter in [50] to even higher powers in the e expansion compared to |46] . The 
only problem is that those masters were computed near four and not six dimensions which is the 
dimension we require them for to complete our master integral determination. This is achieved by 
using the method of |51l [52] . It allows one to relate a d-dimensional integral to a set of integrals 
in (d + 2)-dimensions. The latter set will always be of the same topology as the original lower 
dimensional one but with increased propagator powers. Indeed if the d-dimensional integral has 
P propagators then the {d -|- 2)-dimensional integrals will have {P -|- L) propagators by simple 
power counting where L is the number of loops. In our case if the four dimensional integral 
is a known master then it can be related to the as yet undetermined six dimensional integral. 
Though to do this one has to use the Reduze database to effect the reduction to the unknown 
master and a set of previously evaluated masters. In other words if one builds the system of 
master evaluation from integrals with low numbers of propagators then one can move up the 
tower of unknown masters systematically until all the ones required for the renormalization have 
been found. In the list of masters given in |46j there is enough information in terms of the e 
expansion there in (4 — 2e) dimensions to ascertain all the Feynman graphs contributing to 
the full renormalization of (j2.2p to four loops. To assist with future work we have recorded the 
values of certain integrals in (6 — 2e) dimensions in Appendix A. Rather than present the ones 
used in our efforts, so that others are not restricted to our particular basis of masters, we have 
presented the values for the same topologies given in [36]. These integrals were denoted by Mj 
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and we have given several to quite a high order in e. This was partly because of the spurious pole 
problem. As a simple check on this evaluation of the four loop masters we have followed the same 
procedure at three loops using the three loop four dimensional masters given in [l6] . When these 
were included in the automatic three loop renormalization we correctly reproduced the results 
of [221 [24] . While a simple check we also had to do this exercise anyway as we needed higher 

terms in e at three loops since each three loop integral will be multiplied by counterterms. Hence 
they contribute to the four loop renormalization constants. Finally to effect the renormalization 
within the automatic computation we follow the method of [H]. In this one computes the Green’s 
function in terms of the bare coupling constant. Then after all the graphs have been summed the 
bare variable is rescaled by the coupling constant renormalization constant. This systematically 
introduces the appropriate counterterms automatically and the overall remaining divergence in 
the sum is fixed by the associated unknown renormalization constant. One advantage of this 
approach is that it avoids the subtraction at the level of individual diagrams which is tedious and 
not possible to encapsulate easily in an automatic symbolic manipulation algorithm. Finally, 
with this we have provided all the computational pieces to fully renormalize six dimensional 
theory in the MS scheme at four loops. 


3 Results. 


We are now in a position to formally record the renormalization group functions at four loops 
for various formulations of six dimensional (p^ theory. Throughout all our results will be in the 
MS scheme and we have included an electronic file data with the main results of the article. For 
(12.21) we have 




and 


[-T2 + m]^ + 

o 


- 11 T| + 66 X 2 X 3 - 108Xi - 72 X 5 


i 288 


+ 


-821X^ + 


6078T|r3 - 12564r2r| + 2592C3r2T5 - 9288 X 2 X 3 - 11664^ 


+ 


- 51840C3T3r5 + 61344r3r5 + 20736r7i + 622 O 8 C 3 T 72 - 20736r72] 
-20547T2^ + 1728C3r2^T3 + 185774r|r3 - 31104C3r|r32 - 510960T|r| 


+ 127008C3r22T5 - 23328 ( 4 X^X 5 - 285336r2^r5 + 373248C3r2r33 - 437472r2rJ* 

- 2716416C3T2r3r5 + 559872C4r2r3r5 + 2744064r2r3r5 + 124416C3T2r7i 

- 622080C5r2r7i + 1005696r2r7i + 2457216C3r2r72 - 559872C4r2r72 

- 1005696r2r72 - 1866240C3r3^ + 3005424r3^ + 8315136C3r|r5 - I866240C4TIT5 

- 6763392T|T5 - 5474304C3T3T71 + I492992OC5T3T71 - 7755264T3T71 

- 6096384C3T3r72 + 2239488C4T3T72 + 7755264T3r72 + 1306368C3T| 

- 1321920r| - 7464960C3r9i + 7464960C5r9i - 746496T91 - 7464960C3r98 
+ 746496OC5T98 + 1492992C3r99 - 1866240C5r99 - 3732480C3r92 

+ 373248OC5T92 + 10450944C3r93 - I492992OC5T93 - 2239488C3T94 
+ 746496r94 + 10450944C3r95 - 14929920C5r95 - 746496OC3T96 + 7464960C5T96 

+ 5225472C3T97 - 7464960C5r97] + 0{g^^) (3.1) 




- ^g‘^ +[-11X2+ 24X3] 


432 









+ 

+ 


-821T| + 3222T2T3 - 3060r| + 2592C3T5 - 4536T5 

-I ozzUo 

- 2283 r| + STGCsTlTg + 12428 T|T 3 - 6912C3T2r32 - 14872T2r32 + I4II2C3T2T5 

- 2592C4T2r5 - 20928T2T5 + 41472C3r| - 43392rJ* - 87552C3r3r5 

- 10368C4r3T5 + 122496T3r5 + 13824C3r7i - 69120C5r7i + 59904T71 

+ 34560C3r72 + 31104C4r72 - 59904r72] + 0(5^°) (3.2) 


where Cz is the Riemann zeta function. As a check we have reproduced the three loop results 
of [231IM] and moreover, the non-simple poles in e in the four loop renormalization constants 
correctly emerge. Their residues depend on the poles in e at the previous three loops. A comment 
on our conventions is in order at the outset in relation to other papers. As is usual, [231124], 
we have absorbed the common factor of S{d)/{2 'kY, where S{d) is the surface area of the d- 
dimensional unit sphere, into g^. This factor plays no role in the values of critical exponents. 
Also in comparison with [231 [2Tj our renormalization group functions are defined with an overall 
factor of 2 different. This will be our convention throughout this and later sections. Finally, to 
map our results to [23l [23] and later to isniiii] the sign of needs to be reversed. 

For the mass operator 

O = (3.3) 

the anomalous dimension is 


loig) 


2 I r rp , r>Arp 1 ^^25"^ 

- yff + [-T 2 + 24r3] 


+ 

+ 


-380T| -k 432C3r2T3 -k 71ir2r3 - 864C3r| - 1170T| - 756r5 
-34560C3r2^ + 42635r2^ -k 261792C3r2^r3 - 69984C4r2^T3 -k 364812r|r3 


J 1728 


- 544320C3r2r| -k 419904C4r27f - 1244160C5T2r32 - 200088r2r| 

+ 69984C3r2r5 + 23328 C 4 T 275 - 505872r2T5 + 3825792C3T3^ - 559872C4T| 

- 2488320C5r| - 134784r| - 2975616C3r3r5 - 466560C4r3r5 + 248832 OC 5 T 3 T 5 
+ 2376000r3r5 + 870912C3r7i - 1866240C5r7i -h 1638144T71 + 684288C3r72 

T 0 ® 

+ 559872 C 4 T 72 + 2177280C5r72 - 1638144r72] + 0{g^^) . (3.4) 


The three loop piece agrees with [231124] and we note that we have followed the convention used 
there to include the wave function renormalization constant as part of the operator renormaliza¬ 
tion constant. We have checked that 70 ( 5 ) is equivalent to the /3-function when the expressions 
are reduced to the single coupling theory (|2.1|1 . In the conventions of |23[(24] the relevant relation 
is 

^7?i(5)- 70 ( 5 ) = Pig) (3.5) 

for (|2.in . For completeness we reproduce the independent renormalization group functions for 
(Etd which are 


125 


Pig) = - ^g^ + 5[2592C3 + 6617] 


g 


41472 


[- 4225824C3 -k 349920C4 + 1244160C5 - 3404365] 


g 


+ 0 ( 5 '') 


1492992 


1 
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7 *b) = -75s' + + [ 2592 C,- 51951 -^,^^„^ 
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(3.6) 


0 ® 

+ [lOOSOCs + 18144C4 - 69120C5 + 53449] + 0{g^°) 

which will be relevant to the Lee-Yang edge singularity problem. 

As a final application of the renormalization it is interesting to consider the situation when 
the (j)^ theory is endowed with SU(Nc) symmetry. Such theories were considered at two and three 
loops in [53( [54] as simple models of the strong interactions in four dimensions and denotes 
the number of colours. The motivation was partly due to the six dimensional theory being 
asymptotically free and so it could have parallel properties to QCD. Indeed this was in part the 
starting point for the study in [25] . The other motivation in |53[ IM] rested in the idea that the 
ultraviolet properties of one theory could be regarded as being driven by the infrared behaviour 
of another. So the aim in |53[ 154] was to determine the relevant d-dimensional critical exponents. 
In light of the recent development of the conformal bootstrap method which aims at examining 
the fixed point properties of various scalar theories, we will therefore extend the results of [53[154] 
to four loops. This is in order to provide complementary data for future bootstrap studies. In 
[53] 154] two theories with underlying SU{Nc) symmetry were considered. The first scenario 
was when the tensor is identified as the totally symmetric rank 3 tensor in SU{Nc). In this 
case the group invariants become 

r, = , T3 = 

T 71 = ^[iX,2-8][iV,^-8Ar2 + 256] , + 528] 

Tgi = - 20N^ + 352N^^ - 5120N^ + 26880] 

Tgs = Tge = Tgg = [3N^^ - 16N^^ - 896N^^ + 7296]^ 

Tga = Tgs = Tg-r = [A(,6 + 16A(,^- 1024Ar2^7104]^ 
rg4 = - [iV,6-64Ar4 + 1216iV,2-6784]^ 

Tgg = [5Ar6_72Ar4_ 040^^2 ^3 7) 

where we have made use of the properties of given in |55] when the indices i take values in 
the adjoint representation. In this instance the renormalization group functions are 

3 5 

p(g) = [N^ _ 20] — + [-5iY'^ + 496K^ - 5360] —^ 

i c ^8 Nc ^ ^288iX,2 

-f [211N^ + 62208C3A(,'^ - 27132A(,^ - 20736C3A(,^ -h 1220688iV^2 _ 4395032(3 

-f [-870048(3A)-® -h 1321920(5iV^® - 327893A(,® 14427072(3A(6 559872(4A(6 

- 31570560(5A(,® -b 8142840A(,® - 155416320(3 - 11384064(4 

-b 421770240(5iX,"^ - 112740480iX,'^ -b 1477343232(3iX;^ - 35831808(4iX:^ 

- 1950842880(5 -b 1264882304A(,2 - 7029669888(3 -b 791285760(4 

- 995328000(5 - 5761837824]+ 0{g^^) 

2 4 

l^g) = [-^'c' + 4]4^ + [iVe'-104Ar2+400]^J^ 
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loig) 


+ - lOSeSCsA^;'^ + 17196A^.'^ + 145152C3A^c^ - 296592Ar2 _ 414720^3 

+ [5472C3A^,® - 8640C5iX;^ + 2277N^ - 92160C3A^;® + 15552C4A^.® + 172800C5iVc® 

- 57760A^.® + 804096C3iX;'^ + 393984C4iV/ - 3317760C5A^c^ + 2122784A^.^ 

- 4018176C3A^c^ - 7133184C4iVc^ + 28753920(5- 22816768A^.^ 


+ 7704576(3 + 21233664(4 - 70778880(5 + 60416256] + 0{g^^) 

2 4 

[-N^ + 4]|^ + [lliV,' - IMN^ + 560]^ 

+ [-317N^ + 1728(3A(‘^ + 8664iV^'^ - 27648(3A^^2 _ gogeOiX,^ + 82944(3 


+ 245504] 


9 


1728iV3 


+ [547344(3A(3 - 855360(5A(8 + 362939A^^® - 14211072(3+ 839808(4A(6 
+ 20995200 ( 5 IX,® - 87900801X,® + 166116096(3iX;'^ - 12503808(4iX;‘^ 

- 170449920(5+ 118589664iX^ - 943681536(3iX^ + 64198656(4lX;^ 
+ 423014400(51X2 _ 8129653761X2 1886257152(3 - 110481408(4 


- 89579520(5 +1824079616] 


7464961X'^ 


+ 0 ( 5 '°) 


(3.8) 


where the two loop expressions for f3{g) and 7 ,^( 5 ') are in agreement with [53] when converted to 
the conventions used there. For large IX the /3-function is an alternating series. One interesting 
case is when IX = 4 when the /3-function has a Banks-Zaks type fixed point, CH, in strictly six 
dimensions since 






sr/(4). ^ 

7o [ 9 ) 


-1/ + + |6480C3+2417|4, 

o9 

[-3053376(3 - 58320(4 4786560(5 - 364729]^^^ 0{g^^) 

- - 9 ^ - — 9 ^ + [-48(3-^43] — 

4^ 48 ^ ^256 

(jS 

+ [18864(3 + 10368(4 - 46080(5 + 21907]^ + 0{g^^) 

9216 

3 2 I 9 4 I oroA VI 

-2» + l6« +3[2C3-71.^ 

qS 

[336384(3 55728 ( 4 - 338400(5 -h 335503] 0{g^^) . (3.9) 

y ^ xo 


The second case considered in |S31|S1] was when the symmetry group is SU{3) x SU{3). We 
summarize the relevant background to the construction of that Lagrangian, as it is more involved 
than the previous case, before detailing the extraction of our results. The main ingredient is 
that becomes a complex field and the Lagrangian (12.21) is extended to group to 

L = . (3.10) 

To implement the SU (3) x SU (3) symmetry the coupling tensor is defined by 

(3.11) 
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where the Greek letter subscripts take the values 1, 2 or 3 and eap'y is the Levi-Civita symbol. 
The matrices T* are related to the unit matrix and the SU{3) group generators. In other words 
0 < i < 8 and 


^0 

-^a/3 



(3.12) 


with the remaining eight objects corresponding to the SU{3) group generators. The normaliza¬ 
tion of is chosen, [53l |5l]) so that for the nine objects spanning the group we have 


Tr 




(3.13) 


With this convention we have the simple rule for the full set of T* with 0 < i < 8 that 


^a/3^7A ~ 2^aX^'Yl3 ■ 


(3.14) 


These relations are crucial in evaluating the products of the coupling tensors within the renor¬ 
malization group function determination. We note that unlike [531 M our convention is to 
base the definition of the SU (3) x SU (3) tensor on the SU (3) group generators themselves 
rather than the Gell-Mann matrices used in [53[ [54] . With this construction the fields are in the 
[3,3] -|- [3,3] representation of SU{3) x SU{3). 

To renormalize (I3.10p we constructed the relevant set of Qgraf diagrams to three loops. As 
noted in [54| the number of Feynman diagrams is substantially less than that of the real scalar 
theory (|2.1I) . This is because not all the basic topologies survive when the field is complex. It is 
simple to see this as each line of a graph must have an arrow representing the charge. In addition 
from (|3.10l) the two vertices are represented by the convergence or divergence of three directed 
lines where all lines either have the arrows directed to the point or directed away from the point. 
Therefore, to see which topologies survive from the real theory one merely takes the topology and 
endows its lines with arrows consistent with these rules. If this results in a graph with any vertex 
where not all arrows converge or diverge then it is absent or excluded. It transpires, for instance, 
[54] . that any graph with a subgraph with an odd number of propagators will immediately 
vanish. We have checked this by carrying out the full three loop renormalization explicitly with 
the Lagrangian (|3.10p . Moreover, the renormalization group functions to this order are the same 
as (|3.ip and (j3.2p but with = 0 and Ti 2 = 0. The first invariant corresponds to the basic one 
loop triangle 3-point function. So it is clearly absent, [5l]. Given this observation with regard 
to the renormalization group functions already determined then the four loop expressions will 
be given by excluding those invariants in (13.ip and (|3.2I) which correspond to absent topologies. 
It is straightforward to check that of the set given in ()2.4I1 only Tq 4 passes the test and the 
remaining four loop invariants are excluded from the renormalization group functions. At this 
stage we have only argued for the consequences of the renormalization group functions when one 
has a (j)^ theory of the form (13.101) . For the SU{3) x SU{3) case we need the explicit values of 
the non-zero invariants T 2 , T 5 , T 71 and T 94 . Using the properties for this group discussed above 
we find the values 

T 2 = — , T 5 = 0 , Tji = — 5 T 94 = - (3.15) 

^ 2 32 256 

in our conventions. We note that, [54] . while the two loop nonplanar topology corresponding to 
T 5 passes the arrow test it actually vanishes for the specific SU{3) x SU{3) group. For other 
groups T 5 may be non-zero. These result in 


Pig) 


1 , 11 . 4363 7 

16^ 1152^ 331776^ 
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+ 0(5'') 




io{g) 


+ [-4032C3 - 5760C5 + 10279] 

_ 1 2 _ JL 4 _ 821 6 

24^ 1728^ 497664^ 

+ [1152C3 - 5760C5 +4231]^^^^ + 0 ( 5 '°) 

1 2 1 4 95 6 

- q — - q — - q 

4^ 192^ 3456 

+ [183168C3 - 466560C5 +452171]^^^!^ + 0(g^°) 


(3.16) 


in the MS scheme. 

The two loop expressions given in (|3.16p are in accord with [53t IM] . That they do not agree 
with the three loop results is because (13.101) was renormalized in a different scheme. The scheme 
of [ 54 ] is what is now termed the momentum subtraction scheme which was developed later in 
|56j for QCD. Briefly the coupling constant renormalization constant is defined, |54j . by requiring 
that there are no corrections beyond the tree term when the vertex function is evaluated at the 
fully symmetric point where the values of the square of the external momentum of each leg 
are equivalent. For (I3.inh when the symmetry group is 517(3) x SU{3) the constraint from 
the group theory meant only one three loop graph needed to be determined at this subtraction 
point. Moreover, only the residue of the simple pole was needed for the three loop momentum 
subtraction /3-function. The finite part would be necessary for the four loop correction. That 
7</>(5) agrees at two loop is perhaps somewhat surprising given that only the leading term of 
this is independent of the renormalization scheme. In (j3.inp the agreement at two loop appears 
to be purely coincidental rather as deriving from some property of the two renormalization 
schemes. However, we do actually have a check on the three loop computation which comes 
via the critical exponents computed in |53l [5T| at the Wilson-Fisher fixed point. While the 
renormalization group functions are central to the evolution of Green’s functions for a range of 
scales, the critical exponents at a fixed point are renormalization group invariants. So they have 
the same value in all schemes. In [531 [54] the d-dimensional momentum subtraction scheme 
renormalization group functions were provided. In other words the finite parts of the underlying 
momentum subtraction scheme renormalization constants manifest themselves as 0(e) terms in 
I3{g) and 7 ,/,( 5 ). If one were merely interested in purely six dimensions the 0(e) terms would 
have been absent in the expressions given in |53l 154] . More crucially, if they were excluded 
from a Wilson-Fisher fixed point analysis the derived exponents would not be renormalization 
group invariants. By contrast, in MS there are no 0(e) terms in the coefficients of d-dimensional 
renormalization group functions aside from the 0{g) term of I3{g) which defines the dimension 
of g in d-dimensions. This is because by definition in the MS there are no finite parts in the 
renormalization constants. Therefore, from our MS expressions (j3.16|) we have 

22 4847 6^ 

(jj = 2e - —e'^ --b [-326592C3 - 466560(5 + 422035]^^ + 0(e®) 

2 32 ^4 

V = 3 " + ye" + 16[18C3 + 1]^ + 0{e^) (3.17) 

where there is no O(e^) term in g. Comparing the three loop exponents with the corresponding 
expressions given in [531 [Ml we find exact agreement. This is reassuring since the original 
three loop computation and the present one were in different schemes. However, using the MS 
scheme, where we were able to simplify the computation of the divergence structure of the 3- 
point functions, has allowed us to proceed to a higher loop order for these exponents than would 
be currently possible in a momentum subtraction scheme. 
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4 Lee-Yang edge singularity. 


Equipped with the renormalization group functions we can now study several problems at a new 
order in the perturbative expansion. For the first major application of our results we turn to 
the Lee-Yang edge singularity problem which as elucidated in [29] is related to theory but 
with a purely imaginary coupling constant. In terms of the results in the previous section this 
requires setting = z in the initial Lagrangian to produce the values 

Y 2 = — 1 , T 3 = — 1 , Ts = 1 , Tyj = — 1 , Tgj = 1 . (4.1) 

Alternatively one can use the mapping g ^ ig in (|3.6p . The key quantity of interest for the 
Lee-Yang problem is the critical exponent a which is related to the anomalous dimension of 4> 
of (j2.ip through the hyperscaling law 


[d-2 + r]] 
[d + 2-g] 


(4.2) 


in d-dimensions. With d = 6 — 2e determining the critical coupling constant in d-dimensions 
and expanding the (p held anomalous dimension at that point we hnd 


2 172 e'^ 


+ [- 2783808C3 3779136C4 - 2799360C5 - 3883409] 


4782969 


+ 0{e^) (4.3) 


where the terms to O(e^) are in agreement with [24] but expressed in our conventions for d. 


d 

rj 


a 



3 loop 

4 loop 

3 loop 

4 loop 

5 

- 0.1450 

- 0.1545 

0.3996 

0.3977 

4 

- 0.3173 

- 0.3824 

0.2664 

0.2534 

3 

- 0.4981 

- 0.6805 

0.0913 

0.0562 

2 

- 0.6826 

- 1.0484 

- 0.1458 

- 0.2077 

1 

- 0.8691 

- 1.4860 

- 0.4831 

- 0.5542 


Table 1. Critical exponents rj and a estimates using Fade approximants. 


d 

rj 


a 



3 loop 

4 loop 

3 loop 

4 loop 

5 

- 0.1468 

- 0.1529 

0.3992 

0.3980 

4 

- 0.3280 

- 0.3702 

0.2642 

0.2558 

3 

- 0.5239 

- 0.6446 

0.0862 

0.0630 

2 

- 0.7281 

- 0.9742 

- 0.1540 

- 0.1958 

1 

- 0.9377 

- 1.3583 

- 0.4921 

- 0.5411 


Table 2. Critical exponents rj and a estimates using Pade-Borel method. 


As the Lee-Yang singularity problem stretches across dimensions to d = 1 one has to be 
careful in using the perturbative expansion for large values of e. Therefore, to gain estimates 
for a we have used Fade and Fade-Borel resummation for rj and then evaluated cr through 
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the scaling law. The results for both are given in Tables 1 and 2 where we have used a [2,1] 
approximant at three loops and [3,1] at four loops. We have included the three loop results 
of [21] for comparison. Prom Table 1 it is evident that for dimensions close to 6 convergence 
appears to be present from three to four loops for a. The large discrepancy in rj estimates down 
to d = 4 seems to get washed out in the scaling law. A similar feature is apparent in Table 2 for 
the Pade-Borel application. Though the convergence if anything appears marginally improved. 
The main problem is that the exact values of cr at d = 1 and 2 are not emerging which are — ^ 
and — I respectively. Indeed if anything the four loop estimates in both Tables for these values 
is worse than the three loop ones. This might have been expected as naively setting a value for 
a parameter to be of order 2 in a summed perturbative expansion will mean the larger O(e^) 
term will dominate. One way to handle this is to use a constrained Fade as discussed in [32|. In 
that method the two exact values for a are included in the derivation of the rational polynomial 
of the Fade approximant. We have carried this out for the four loop estimate of a which is given 
by 


1 1 79 

a = -e-+ [15552C3 - 104451- 

2 6 972 ^ 057464 


+ [-2503872C3 + 3779136C4 - 2799360C5 - 4047533] 


25509168 


+ 0{€^) . (4.4) 


The results are given in Table 3 where we have reproduced the constrained three loop [3, 2] 
results of [32] and given our four loop [4,2] Fade estimates. The constraints are included 
for completeness. Results from other methods are included for comparison. These include a 
strong coupling expansion |32|; as well as two Monte Carlo methods which are based on critical 
behaviour in problems seemingly unrelated to the Lee-Yang singularity problem. These are 
termed (lattice) animals and fluids with the former originating in polymers in a solvent and the 
latter related to pressure in fluids where there is a repulsive core, [solEI]. The final column in 
Table 3 are recent results from a conformal bootstrap analysis, [33|. With the inclusion of the 
exact results for low dimension in the Fade approximant not only is there better convergence 
for d = 3 from three to four loops but there is remarkable agreement with the values from |32j . 
For the other methods the four loop estimates lie within error bars except compared to the 
d = 3 value for the fluids method. In light of this it is worth commenting on why we presented 
Tables 1 and 2 in the hrst place. This is partly to make contact with [23|. More crucially given 
that we will be using Fade approximants later for other problems down to low dimensions it 
is important to be aware of the potential limitations of the technique in the absence of known 
exact two dimensional conformal field theory exponents. As exact results are known for d = 1 
and 2 here we can gauge how far off estimates may be for these dimensions. The d = 3 values 
for a in Tables 1 and 2 are perhaps not reliable but those for d = 4 and 5 are in keeping with 
those of the other methods listed in Table 3. 


d 

3 loop 

4 loop 

Ref [32] 

Ref [30] 

Ref [31] 

Ref [33] 

5 

0.3989 

0.3981 

0.401(9) 

0.402(5) 

0.40(2) 

0.4105(5) 

4 

0.2616 

0.2584 

0.258(5) 

0.2648(15) 

0.261(12) 

0.2685(1) 

3 

0.0785 

0.0747 

0.076(2) 

0.0877(25) 

0.080(7) 

0.085(1) 

2 

- 0.1667 

- 0.1667 

- 0.166(5) 

- 0.161(8) 

- 0.165(6) 

- 0.1664(5) 

1 

- 0.5000 

- 0.5000 

- 

- 

- 

- 


Table 3. Critical exponent a estimates using constrained Fade approximant and comparison 

with [30l|3ll[32ll33|. 
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5 Percolation. 


We now turn to the application of the renormalization to the percolation problem which requires 
the evaluation of the tensors (j2.4p for a specific configuration. The percolation problem is 
described by the replica limit in the {N + l)-state Potts model [34j which in the case of (12.2p 
corresponds to a special designation of the coupling tensor A straightforward way of 

mapping to the Potts model was provided in ISZ! which involves a set of vectors, e^. These (A^+1) 
vectors describe the vertices of an Wdimensional tetrahedron and allow one to decompose 


as 


Jijk ^ 


N+l 


V e® 


~'(y. OL OL 


OL=l 


In order to represent the tetrahedron the vectors must satisfy the following relations, 

N+l N+l 

T 4 = 0 . T = (AT + 1)0 

a=l a=l 

for sums over the (N + l)-dimensional label and 


N 

E 

i=l 


~ (-^ + l)<5a/3 — 1 


(5.1) 


(5.2) 


(5.3) 


for summations over the original indices denoted by i. It is the form of the final relation 
which means that the evaluation of the underlying tensors (12.41) requires special care. In |24| a 
diagrammatic method was outlined to handle the lower rank tensors. However, we have written 
a Form routine to reproduce the evaluations given in [21] and then applied it to the cases Tg^ 
for 1 < m < 9. Such a systematic path seems more appropriate since the diagrammatic method 
is tedious at three loops as it involves seven sums but manageable for only two tensors. At four 
loops the nine summations for nine independent tensors is not straightforward. For arbitrary N 
we have 


T2 = 

= {N - 

f l)^[iV 

-1] , Ta = (iV + l)2[iV-2] 



n -- 

= m 

+ lf- 

6(IV + l) + 10](iV + l)^ 




Tn -- 

-- 

+ 1)^ - 

9(A^ + l)2 + 29(iV + l)- 

32](iV + 1)® 



T 72 = 

= m 

+ lf- 

6(A^ + l) + ll][iV-2](IV 

+ 1)® 



Tgi -- 

-- 

+ 1)^ - 

9(IV + 1)2 + 30(iV + 1)- 

35][iV-2](iV 

+ lf 


Tg2 = 

= Tge 

= Tgs 

= [{N + if - 9{N + if 

+ 30(iV + 1) 

-38] [AT 

-2](iV + l)® 

Tg3 = 

= 2g5 

= Tgj 

= [{N + if - 9{N + if 

+ 30(iV + 1) 

-37][A^ 

-2]{N + lf 

Tg4 = 

= [{N 

+ 1)"- 

12{N + if + 57{N + if 

- 125(A^ + 1) 

+ 106] (IV + if 

Tgg = 

= 

+ lf- 

5(IV + l) + 10][iV-2][Ar 

-3]iN + lf 




(5.4) 

where we present the expressions in the same format as |24] and have included the known values 
from pi] for completeness but with T 72 factorized further. Although there are two instances 
of three tensors giving the same value for this tetrahedron conhguration, it is clear from the 
underlying Feynman diagram dehning the tensors that the graphs themselves are topologically 
distinct. This should not be a surprise since in QCD, for example, when one examines high 
loop diagrams different topologies can have the same combination of colour group Casimirs 
multiplying them. Taking the N ^ 0 replica limit gives the values we require for the evaluation 
of the critical exponents for the percolation problem. We have 


T 2 — — 1 , Ts — — 2 , T 5 — 5 , T 71 — — II , T 72 — — 12 

Tgi = 26 , rg 2 = 32 , Tga = 30 , rg 4 = 27 , Tgg = 36 


(5.5) 


16 


for the independent tensors. 

With these values we have computed various critical exponent in powers of e to 0{e^). Using 


7 


671 


Pig) = -off - off - 7^9 + 


+ 

1 


8 

121109 


288 ^ 


414031 93 ■ 

41472 ~ 16*’^ 


ff 


651 595 84156383 

1728 ~ ~ 1492992 


ff" + 0{g^^) 


/ N - 2 37 4 

7</.(ff) - ^2^ 432^ 


29297 5 ■ 

62208 ~ 


70 (ff) = 


+ 

^ 2 
2 ^ 

+ 


225455 233 33 55 ' 

82944 ~ 


ff" + 0 (ff'") 


47 4 ^ 
48^ + 


3709 3 ^ 

-h -C 3 

864 2^"^ 


ff 


18486131 20027 33 15 

746496 864 ~ 


ff" + 0 (ff'°) 


(5.6) 


they are 


r? =-e — 

' 21 


824 
9261' 


+ 4[290304C3 - 93619] 


4084101 


+ 2[286336512C3 + 384072192C4 - 1493614080C5 - 103309103] 


1801088541 


+ 0 (e") 


go = 


7 = 


u = 


4 

- - 


710 


+ [925344C3 - 235495]- 


3087 ' ^ ^1361367 

+ [603983520C3 +1224230112C4 - 5334336000C5 - 157609181] 


1200725694 




565 


+ [-925344C3 +408997]; 


3087 ' " '2722734 

+ [-933950304C3 - 12242301 12 C 4 + 5334336000C5 + 302378687] 
1 5 589 


2401451388 


+ 0 ( 6 ^) 


+ 0{e^) 


+ 


+ + 


+ [-1614816C3+ 716519]- 


w = 2e — 


2 ' 42 9261 ' '16336404 

+ [344397667 - 1344827232C3 - 2136401568C4 + 10028551680(5] 
1342 


14408708328 


+ 0{e>) 


-e^ + [62496(3 + 40639]; 


441 ' '7203 

+ [248046624(4 - 702654624(3 - 1209116160(5 - 317288185]- 


, + 0(e®) .(5.7) 

'19059138 y J \ J 

The exponents rj, rjQ and u are obtained from the corresponding renormalization group function 

and 7 and u are deduced from the scaling relations 

-1 


go = i' - 2 + r] 


7 = (2-ff)i^ 


(5.8) 


With these we have repeated the exercise of the previous section to obtain Fade and Pad+Borel 
estimates of various critical exponents following the method of [23] . There estimates were found 
for rj and ryo and then values for 7 and u were obtained from the scaling laws. Our results 
are contained in Tables 3 to 7 with those at three loops agreeing with [23j and included for 
comparison with the four loop estimates. The three loop estimates for u are in accord with 
those given in 


Overall a similar feature emerges as for the Lee-Yang edge singularity estimates in that 
down to four dimensions there is reasonable convergence but below this the results are not as 
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d 

rj 

Vo 

7 

u 

OJ 

5 

- 0.0569 

- 0.3097 

1.1773 

0.5723 

0.7910 

4 

- 0.1186 

- 0.6319 

1.4250 

0.6726 

1.5155 

3 

- 0.1812 

- 0.9566 

1.7812 

0.8166 

2.2326 

2 

- 0.2442 

- 1.2822 

2.3328 

1.0395 

2.9473 


Table 4. Critical exponent estimates using Fade approximants to three loop expressions. 


d 

V 

Vo 

7 

V 


5 

- 0.0594 

- 0.3192 

1.1834 

0.5746 

0.7085 

4 

- 0.1338 

- 0.6885 

1.4764 

0.6919 

1.1590 

3 

- 0.2215 

- 1.1048 

1.9893 

0.8955 

1.4775 

2 

- 0.3222 

- 1.5678 

3.0782 

1.3256 

1.7150 


Table 5. Critical exponent estimates using Fade approximants to four loop expressions. 


reliable. This situation was ameliorated by exploiting known results in two dimensions and 
using this as a constraint or boundary condition on the Fade approximant. Therefore, we have 
followed this procedure again and constructed constrained Fade approximants to 7 and a; 
from their respective exact two dimensional values of H and 2, [59] . The results of this 
exercise are given in Table 8 . There the two dimensional values of rj and the exponent j3 agree 
with their known exact values. The latter exponent as well as a and r, which are deduced from 
hyperscaling laws from the previous columns in the table, are included for comparison with 
results from other methods. In this respect Table I of [60] gives a comprehensive summary of 
estimates for these exponents. To compare we have included the results of the computation 
of v, 7 , rj and /3, [60], in Table 9 which used a high temperature series method. Though it is 
worth noting that a more recent study, [61], has obtained estimates for 7 which are equivalent 
to the those of the high temperature series of [60] . In m the series was extended beyond the 
15th order of [60]. Examining the results given in Tables 8 and 9 on the whole the constrained 
Fade estimates are in reasonable consistency with the central values of [60] . Ferhaps more 
significantly the three dimensional estimates from the perturbative approach are in line not 
only with [60] but results from other methods as is evident from Table I of [60]. There the 
perturbative estimate of 0.34(4) was quoted for /3, for example, but that of Table 8 is more in 
line with other methods now. One exponent not covered by the summary table of [60] is the 
correction to scaling exponent to. Two studies of iv using Monte Carlo methods are given in 
[62] and [58] which give results in three and four dimensions respectively. These are 1.62(13) 
and 1.13(10). Other estimates in three dimensions are 1.61(5), [63], and 1.77(13), [64] . Our 
estimates in Table 8 are in remarkable agreement in three dimensions and within the error of 
[58] in four dimensions. With the results in Tables 8 and 9 and the close tally it gives support to 
the earlier observation that the usual Fade approximant is only really competitive down to four 
dimensions. Below that the summation fluctuates as is apparent from the four loop estimates 
and does not capture the exact two dimensional picture. Hence results down to four dimensions 
should only be considered. For the remaining two exponents a and r we note that the estimates 
for r are in good agreement over all dimensions in comparison with Monte Carlo estimates given 
in Table I of m- An estimate for a of around 0.452-0.454 is given there too but only in three 
dimensions. So our value is below the central value. The remaining two columns in Table 9 
correspond to estimates for a and r from [63l [65l [ 66 ] . Again the constrained Fade estimates are 
not dissimilar to the central values. 
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d 

77 

Vo 

7 

u 

OJ 

5 

- 0.0578 

- 0.3122 

1.1788 

0.5729 

0.7540 

4 

- 0.1229 

- 0.6431 

1.4346 

0.6758 

1.3880 

3 

- 0.1903 

- 0.9800 

1.8097 

0.8262 

1.9937 

2 

- 0.2588 

- 1.3199 

2.4057 

1.0651 

2.5870 


Table 6. Critical exponent estimates using Pade-Borel method for three loop expressions. 


d 

V 

Vo 

7 

V 


5 

- 0.0582 

- 0.3160 

1.1814 

0.5740 

0.7539 

4 

- 0.1253 

- 0.6631 

1.4517 

0.6831 

1.3966 

3 

- 0.1968 

- 1.0464 

1.9096 

0.8693 

2.0170 

2 

- 0.2717 

- 1.4503 

2.7656 

1.2174 

2.6320 


Table 7. Critical exponent estimates using Pade-Borel method for four loop expressions. 


We close this section by recording the renormalization group functions for the N = 2 case. 
While this is not directly related to the percolation problem it does correspond to a specific 
Potts model and so we give the relevant expressions for completeness here. The main reason for 
treating this case specially lies in the nature of the group invariants T„ given in (|5.4I) . From the 
explicit values most vanish at = 2. If one analyses the underlying Feynman graph which each 
invariant relates to, then the non-zero invariants have the same feature as the SU{3) x SU{3) 
theory analysed earlier. Though they have no other connection aside from the graphical one. 
This similarity is that when is non-zero the corresponding Feynman graph has no subgraph 
with an odd number of Feynman propagators. If any Tj = 0 then there is at least one subgraph 
with an odd number of propagators. This is clearly the case for Ts which is the one loop triangle. 
So in (j3.ip a large number of terms are absent. Moreover, the non-zero values of are 


Ta = 9 , Ts = 81 , Tn = 729 , = 6561 


(5.9) 


which are all powers of 3 in contrast to ()3.10p where the non-zero invariants were powers of 1/2 
and Ts is non-zero here. With these values the renormalization group functions are 


/ 3 ( 5 ) = - ^9 - - 


747 


9 


+ 9[2592C3 + 106271-^^2 


-F 243[-25248C3 - 864(4 - 23040(5 4607] 


2048 


+ 0{g^^) 


33 


= --9- ^9^ + 3[2592(3 - 5357] 


256 


-F 243(3104(3 - 288(4 - 7680(5 + 4077] 


1024 


+ 0(5'°) 


92 27 4 1917 

70(5) = - 2^ - r ^9 - 


-h 9(906336(3 + 23328(4 - 1866240(5 -h 1174907] 


1024 


+ 0{g^^) (5.10) 


which can be used for a Wilson-Fisher fixed point analysis. 
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6 0{N) symmetric theory. 


We now turn to another version of (j)^ theory which is the one endowed with an 0{N) symmetry. 
It has been considered recently in [THl |20l [26l HTJ [28] in the context of the conformal bootstrap 
programme as a way of accessing five dimensional quantum field theories with a conformal 
symmetry. The basic Lagrangian in this case is 

L ^ ^ + ycr^ ( 6 . 1 ) 


where there is an 0{N) multiplet of fields together with a single scalar field, a. To ensure these 
fields produce a renormalizable Lagrangian in six dimensions there are two massless coupling 
constants gi and g 2 - For ease of comparison with (TS] [20] we use the same notation as those 
articles. As the main interest in this section is the extension of the three loop analysis to the next 
loop order and as the computational techniques have already been described we will mention 
only those features which are new to this calculation. First, with the extra fields we have carried 
out the full computation rather than identify the group structures of section 3 with those of the 
0{N) symmetric Lagrangian as was noted in [20]. While this may seem to be inefficient it is 
actually a necessary first step in the derivation of the renormalization of the mass operators of 
the fields in (|6.1I) which we will discuss later. What is worth noting is the number of 2-point 
graphs which are required for the renormalization of (16.11) . For there are 1, 5, 48 and 637 
one, two, three and four loop graphs respectively. The corresponding numbers for the a 2-point 
function are 2, 7, 60 and 723. The resultant four loop field anomalous dimensions are 


j<p{gi,g2) = 


-^ + 
6 


EL 

J 432 


+ 


+ 


- llNgl + 26gl+A8gig2-ng^2 

UN'^gf - 232gfN + 5184C35i - 9064c/f -h 26i6Ngfg2 - 3264c/? c /2 

- 386Ngigj + 5184C35?52 “ 'i-'^762gfgj + M2gigl + 327(7|] 

1296C3iV^5i + 3N‘^gl + 46656C3A^^c/i + 2Ul2N'^gl + 3649536C3A^c/i 

-F 1026432C4A^5i - 3732480C5AIc/? - 1600648iV5? - 1275264C35? 

-F 1306368C45? - 7464960(55? + 90959445? - 15552C3iV^5?52 
-F 52452iV2^?52 - 2799360C3iV5?52 - 839808(4A^5?52 + 3MM32Nglg2 
+ 995328(35?52 + 3359232(45?52 - 27842405?52 + 1296(3iV25f5i 

- 387m^g\gl + 1034208(3iV5f52 “ 116640(4iV5?5| - 28M3lQNgigl 

+ 3037824(35^51 + 233280(45^51 - 14929920(55t5l + 139290645? 5 ! 

+ 77760(3iV5?5i + 3bbUNglgl - 1772928(35?5i + 2239488(45?5l 
+ 19104965? 5 I - 1296(3iV5?5l + 409511V5?5? + 1648512(35?5? 

+ 886464(45?5l - 3732480(55151 + 10566205? 5 I - 342144(35i5l 


d 

V 

7 

V 


iv 

CJ 

T 

5 

0.5746 

1.1817 

- 0.0565 

0.8457 

0.7178 

0.4933 

2.4171 

4 

0.6920 

1.4500 

- 0.0954 

0.6590 

1.2198 

0.4742 

2.3124 

3 

0.8968 

1.8357 

- 0.0470 

0.4273 

1.6334 

0.4419 

2.1888 

2 

1.3333 

2.3888 

0.2083 

0.1389 

2.0000 

0.3956 

2.0549 


Table 8. Critical exponent estimates using constrained Fade approximant of four loop results. 
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d 

V 

7 

5 

/? 

a 

r 

5 

0.571(3) 

1.185(5) 

- 0.075(20) 

0.845(5) 

- 

2.412(4) 

4 

0.678(50) 

1.435(15) 

- 0.12(4) 

0.639(20) 

— 

2.313(2) 

3 

0.872(70) 

1.805(20) 

- 0.07(5) 

0.405(25) 

0.445(10) 

2.190(2) 

2 

1.3333 

2.3888 

0.2083 

0.1389 

- 

- 


Table 9. Results for critical exponents from [HO l [63 | [65 } [ 66 ], 


and 


- 279936(45152 + 4596125i5| + 68688(353 + 23328(452 


- 204484^2 


5? 


+ 0(5^) 


J 6718464 


( 6 . 2 ) 


1<t{9i-,92) — “ -^5i+52 TX+ 2A^5f + 48iV5f5r2 — lliVg'^(72 + I 35 I 


432 


+ 


+ 


- 27Q2N‘^g\ + blSiCsNgf - 8560Ngf + 1152N‘^gfg2 + W56Nglg2 

+ 31V^5i 52 + 12960(3^75^52 - 26646iV5i52 - 1560iV5f5^ + 952Ng‘fg^ 

+ 25926!,!-51955.!]^ 

- 41472(3iV^5? + 54266A^^5^ + 1897344(3 513216C4N^gf 

- 1866240(5N'^gf - 605816iV25? - 238464(3iV5® + 653184(4iV5® 

- 3732480(5iV5? + 3883280A^5? - 8064Ar3^''52 - 2488320(3iV^5('52 

+ 4171512Ar2c,Jc,2 + 1679616(4A^5[52 - 1008768A^5j'52 + 2b02(^N^glgl 

- 364N^glgl + 1642240(5N‘^g\gl - 233280(^N‘^g^gl - 2324662N‘^g^gl 
+ 3535488(3 iV5?5^ + 746496(4 iV5?5i - 14929920(5 A^5)*5i 

+ 10883728A^5)*5| - 233280(’iN‘^g\gl + 416016N‘^g\gl - l026432('iNg\gl 
+ 2709360(^Ng\gl - 240816Nglgl + 12060(5N‘^g\gt - lOlOlN^gigt 

+ 1031616(3iV5^52 + 1283040(4iV5i52 “ 7464960(5A^5t52 

+ 6462626iV5^5| - 108864(3A^c/?5| + 289416A^5?5i + 196992(3iV5f5f 

- 03312(iNglgl - 306628Nglgl + 272160(35^ + 489888(45^ 


- 1866240(55^ + I 4431235 I 


1 


6718464 


+ 0 {gf) 


(6.3) 


where the order symbol refers to any combination of the two couplings. As we do not use the 
method of subtractions but follow the automatic renormalization algorithm of m the cancella¬ 
tion of the double, triple and quadruple poles in e act as a useful computational check. Moreover, 
we correctly reproduced the three loop results given in m- 

For the two /3-functions, /3i(5i,52) and / 32 ( 5 i, 52 ), we again do not need to compute any 
nullified 3-point vertex graphs but instead use the method outlined in section 2 where the 
propagators of each field had an extra term. In terms of graphs this means we avoid computing 
6455 graphs at four loops for gi and 6681 for g 2 . However, in expanding the propagator within 
a graph to reproduce the corresponding nullified 3-point vertices, one has to label the extra 
term of the propagator with the coupling constant 51 . This is because this is the only vertex 
with two fields. Equally the additional term for the a propagator is labelled with g 2 . Once 
this step has been achieved the renormalization process outlined for the basic Lagrangian is 
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followed. The result is the two /^-functions 


/3i(5'i,5'2) = - N gl + Sgl + I2gig2 - gl 

+ [- mNg\ - 5365^ + l^2Nglg2 - 3605?52 - UNglgl - Q2Sglgl 

- 2%! (/I + 135 ! 

+ [3662iV2^f 129600C3 A^c/? - dOhSSiV^^ + 20736(35? + 2513605? 

- 36iV2g?52 - 155520(3 A^5?52 + 124704iV5?52 + 186624(35?52 + 180005?52 
+ ^N‘^g\gl + 12960(3lV5?5i - ^^^^^Ng\gl - 41472(35?5i + 3584805?5i 

- 4560iV5?5i + 124416(35?5i + 977765?5^ + ^^^Nglg^ + 62208(35?5| 

+ 99605?52 - 31104(3515^ + 336125i5? + 2592(35? - 51955?] 

+ [93312(3iV^5? - 12310Ar^5? + 21959424(3^^5? - 5365440(4iV^5? 

- 1866240(5 IV^^? - 115353841V2^? - 172969344(3Ar5? + 14183424(4 IV 5 ? 

+ 111974400(5 iV 5 ? + 12401088iV5? + 31290624(35? + 1492992(45? 

- 82114560(55? - 1046803845? - 31104(3iV35?52 + 42481V^5?52 

- 17262720(3 iV25?52 + 2799360(4iV25?52 + 119981521V2^?52 

- 2967321 6(3 iV 5?52 + 3919104(4lV5?52 + 89579520 ( 51 V 5?52 

- 93820080iV5?52 - 89268480(3 5^52 + 12877056(45?52 - 78382080(5 5^52 
+ 59029445?52 + 2592C3N^gfgj - 35AN^gfgj + 3403296(3N^gfgl 

- 233280(4 iV^ 5 ? 5 i - 4985248iV25?5? + 27454464(3 iV 5 ? 5 | 

- 12223872(4 iV 5 ? 5 | - 59719680(51V5?5? + 25092064iV5?5i 

- 176380416(35?5i + 18382464(45?5i + 82114560(55?5i 

- 1096781925? 5 I - 342144(3iV25?5? + 5766481V25?5? - 17324928(3lV5?5? 

+ 839808(4 iV 5 ? 5 ? + 111974400(5 iV 5 ? 5 ? - 41588256iV5?5? 

- 70232832(35?5| + 559872(45?5? + 44789760(55?5? - 1462151525? 5 ? 

+ 12960(3iV25?5? - 191011V2^?5? + 3478464(3iV5?5? - 3195936(4iV5?5? 

- 7464960(5 1V5?5? + 7579786iV5?5? - 85784832(35?5? + 5412096(45?5? 

+ 3732480(55?5? - 340306885? 5 ? - 2255040(3iV5?5? + 1119744(4iV5?5? 

+ 1918896Ar5?5? - 24292224(35?5? + 8957952(45?5? - 55987200(55?5? 

+ 170966165? 5 ? + 196992(3iV5?5? - 93312(4lV5?5? - 3065281V5?5? 

- 6371136(35?5? + 4665600(45?5? - 174261445? 5 ? - 4494528(35x5? 

- 4199040(4515? + 22394880 ( 5515 ? - 99443045i5? + 272160(35? 

+ 489888(45? - 1866240(55? + 14431235?] + 0 ( 5 ?') (6.4) 

and 

^ 2 ( 51 , 52 ) = - ^ + 4iV5? - ^5 i52 + 35 ? ^ 

+ — 24A^ 5 ? — 322N 5?52 — 60N 5 ? 5 ? + 31A^5?5? — 1255? 

+ [27696Ar2^? + 34224Ngl - 38474N‘^gfg2 + 5184(3Ar5?52 + 59408iV5?52 
+ ll304N^gfgl + 62208(3iV5?5i + 2h296Ng\gl - 789N^g\gl 
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+ U064C3Ngfgi + 127890iV5f53 _ 20736C3Ngfgl - 8688Ngfgl 
- m72Nglgl + I296OC352" + 33085(72"] 4^^ 


+ 


1088640C3A^^5? - 1031208iV35? - 8771328C3iV^5? - 3359232C^N‘^gf 
+ 6915984A^25? - 6117120C3A^5? - 559872C4A^5? + 11197440(5iV^? 

- 20404128iV5? - 1021248C3lV^5i52 + 7064781V^5?52 + 3856896(3 
+ 8071488(4lV^5?52 - 13063680(5iV2^?^2 - 26286776iV25?52 

- 34763904(3iV5?52 - 4385664(4A^5f52 + 1866240(5A^5?52 

- 401008iV5?52 + 279936C3N^glgl - imSAN^glg^ 

- 12192768C3N'^glg2 + 1399680(4+ 22394880(5 
+ b773Q32N‘^glgl - 39688704(3iVc/Jc/l + 8398080(4iV5l52 
+ 44789760(5iVc/Jc/l - 67219056A^5[5i - 23328C,3N^g\gl 
+ 999m^glgl + 4388256(3iV"(7?5i “ 3172608(4A^2(7?(7i 

+ 10267192iV25^5f - 56619648(3iVc/?c/^ + 7744896(4A^5i5i 
+ 7464960(5lV(7?(7i - 9887792Nglgl - 1477440(3 

+ 559872(4- 3730536iV^5f5| - 8398080(3iVc/f+ 559872(4iVc/^c/^ 

- 55987200(5iV5f52 " 8554440iV(7?5| + 44064(3Ar25f(7f + 2248171V2^t(7i 

- 22187520(3iV5^c/f - 956448(4lV5f5i + 14929920(5 

_ 44490442iV5^5f - 9782208(3+ 1959552(4Ar5?c/f 
+ 22394880(5 A^5 ?c/^ + 3707040A^5?5^ + 1223424(3A^5f c/^ 

513216(4A^5i 5| + 1351296Ar525^ _ 12677472(35^ + 1049760(45f 

1 


+ 3732480(55^ - 102130955^ 


J 4478976 


+ 0 ( 5 ,"") 


(6.5) 


As with the two field anomalous dimensions we have reproduced the three loop expressions given 
in m and equally recovered the expressions given earlier for (|2.ip when N = 0. This comparison 
is made with reference to the comments on our conventions. In particular the leading term of 
each /3-function derives from our choice of (i = 6 — 2e and the factor of 2 in our definition of the 
renormalization group functions in comparison to |2n( [2T] . 

To complete the four loop renormalization of (16.ip we compute the renormalization group 
functions associated with the two mass operators which we will denote by 


Oi 




( 6 . 6 ) 


These two operators have the same canonical dimension of 2 and therefore mix under renormal¬ 
ization. Here we note that our canonical dimension convention derives from the dimensionality 
of the associated coupling constant of the operators in a Lagrangian. This is an important 
distinction and is motivated in part for later discussion in relation to checks with the large N 
expansion of critical exponents. In other words denoting the bare operators with a subscript o 


0,0 = Z^jOj (6.7) 

where Zij is the mixing matrix of renormalization constants. These produce a mixing matrix of 
mass anomalous dimensions denoted by 7,j(5i, 52 )- To extract the renormalization constants one 
ordinarily inserts the operators into separate (/>* and a 2-point functions of all possible 1-particle 
irreducible Feynman graphs and follows the normal procedure. However, as the operators do not 
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involve derivatives there is no complication with mixing into total derivative operators. More¬ 
over, this means that for the determination of the MS renormalization constants the operators 
are inserted with no momentum flowing in or out of the vertex itself. This is a standard method 
and reduces the problem to a simple 2-point function computation. In the context of the un¬ 
derlying (j)^ interaction this leads to a similar computational simplification which we exploited 
before. Again one need not generate any more Feynman graphs than those already used for 
the wave function renormalization. Similar to the coupling constant we expand each (/>* and a 
propagator as if there was a respective mass present using the mapping 

1 1 

¥ ^ ^ (A:2)2 • 

The additional complication here is that one has to to label the 0{m?‘) term to indicate whether 
that insertion is from a (/>* or cj field mass term. By contrast if one instead evaluated the 3- 
point functions with a nullified operator insertion then there would be 4 one loop, 38 two loop, 
722 three loop and 13136 four loop graphs to determine for the mixing matrix in total at each 
loop order. Following the procedure for renormalizing a mixing matrix we find the elements of 
to four loops are 


711 ( 51 , 52 ) = % + [-22iV5i - 134c/f - 2,^9192 + 652 ] 7 ^ 

6 L -I ZLb 

+ [803iv25t -h 15552C3A^5t - 4016Alc/f -h 2592C35t -h 31420c/f - 7776C3A^5?52 
+ 225mgfg2 + 15552C35?52 - 296453^2 + 3926Nglgl - 5184(35?51 
+ 18512gfgl - 2859gigl - 51gl] 

-F [l6848C3iV35? - 8322N^gf + 2507760(3A^^5? - 734832(4 

- 1366196iV2c/? - 21591360(3-h 1691280(4iV5? -F 14463360(5 
-F 1064726Ngf + 3911328(35? -F 186624(45? - 10264320(55? 

- 130850485? - 909792(31V2^?52 -b 139968(4fV^5?52 + 370827iV25?52 

- 225504(3 iV5?52 -b 1959552(41V5?52 -b 2799360(51V5?52 

- 73120561V5?52 - 8040384(35?52 + 979776(45?52 - 8398080(55?52 
+ 39629525?52 + 104328(3Ar25?52 _ 31177N^gfgi - imiOAC^Ngjgi 

- 991440(4iV5?5l - 2799360(5Ngfgj + 377U32Ngtgl - 9916992(35t5l 

+ 1877904(45?5l + 8864640(55t5l - 117671425?5i - 524880(3iV5?5l 

- 69984(4lV5?5l + 2799360(5^5?5l - 2267862Ngfgl + 1741824(35?5i 

- 979776(45?5l “ 5598720(55?5i + 326645? 5 I + 76464(3lV5?5l 

- 104976(4iV5?5l - 3032991V5?5l - 3335904(35?5l - 443232(45?5l 
+ 7464960(55?5l - 24287085? 5 I - 73872(35i5l + 139968(45i5l 

1 

+ I 3445 I 515 I - 11016(351 - 11664(451 + 3659651 

J 1679616 

712 ( 51 , 52 ) = + - 25 ? - I 85152 - 35 I 

+ [i 154A^5? + 14265? - 992iV5?52 + 18225?52 + 141iV5?5i + 864(35?5l 

+ 14305?5l + 864(35151 + I 42 O 515 I - 2 I 5 I] ^ 

-b [45360(3A1^5? - 42967iV^5? - 365472(3iV5? - 139968(4^5? + 288166iV5? 


24 



i2i{gi,92) — 


and 

722 ( 5 '!, 52 ) = 


+ 


- 254880C351 - 23328C45? + 466560C55? - 850172^^ - 27216C3A^^5f52 
+ 18117N‘^glg2 + llU56C3Ngfg2 + 233280C4A^5?52 - 311040C5iV5f52 

- 856464iV(/f 52 - 1161216C35f52 - 186624C45f52 + 155520C55f52 
+ 881285^52 + 3888C3A^'5^5i - WSbN^gjgl - 260928C3Ngjgi 

+ 66096C4Ngtgl + SimoaNgfg^ + 130289Ngtgl - 1468800C35t5i 
+ 124416C45t5i + 2488320C55^5i - 23602545^^2 ^ iii888C3iV5353 

- 42768C4iV5353 + 135612^^5? 5 ^ - 506304C35?5f + 62208C45?5i 

- 777600C55?5i + 1074185? 5 ! - 7776C3lV5?5l - 19521iV52^4 

- 251856C35?52 + 31104C45?5l - 466560C55?52 " 2094575 ?52 

- 247104C3515! + 50544C45i5i - 363531515^ - 37584C352 + 11664C452 

+ 511655^^ 


i 186624 


+ 0(5?°) 


9j 

2 

+ 


+ 


7 Ngl - 205? - 545152 - 2 gl 


9j 

72 


-99N‘^gf - 7776C3Ngf + 8798A^5f + 5184C35f + 34765? - 4896iV5?52 
+ 175325?52 - 250Ngjgl + 10368C35?5i + 100545?5i + 5184C35i5i 

+ 864515? - 2592C35I + 28015?] 

J 10368 

-1296C3iV°5? + 513A^°5? - 312336(3A^^5? + 69984C4iV^5? + 202501iV2^6 

- 1161216C3iV5? - 559872C4A^5? + 2799360C5A^5? - 14297861V5? 

- 1039392(35? + 209952(45? - 466560(55? - 10750285? 

+ 42768(3 1V2^5^2 - 1004701V2^5^2 + 1492992(3iV5?52 - 209952C^Nglg2 

- 933120(5iV5?52 - 665130iV5?52 - 4043520(35?52 + 139968(45?52 
+ 466560(55?52 - 6475445?52 - 4536(3iV25?52 + m93N^g\gl 

- 504144(3iV5?5i - 58320(4iV5?5i + 3732480(5iV5?5i - 966856Ngfgl 

- 3506976(35?5i + 349920(45?5? + 3732480(55?5i - 61031865?5i 
+ 76464(3iV5?5i - 151632(4iV5?5i + U7648Ngfgl - 2462400(35?5i 
+ 373248(45?5? - 2332800(55?5i - 6083765?5? - 89424(3iV5?5? 

+ 46656(4 1V5?5? + 67895iV5?5? - 987552(35?5? + 326592(45?5? 

- 2332800(55?5? + 6675425?5? - 261792(3515? + 198288(4515? 

- 738288515? - 187272(35? - 174960(45? + 933120(55? 

-4143465?! + 0(5?°) (6.9) 


J 559872 


—Ngf + 55 ? 


1 

J 12 


1 

J 216 


+ 


-80iV5? - 301V5?52 + 26Ar5?5i - 975?J 
-20618iV25? + 5184(3Af5? + 27800iV5? + 11304Ar25?52 + 62208(3iV5?52 

- 14064iV5?52 - 11851V2^?5i + 28512(3lV5?5i + 1540381V5?5i 

- 31104(3lV5?5? - 11496iV5?5? - 98841V5?5? + 18144(35? 

+ 522255 ? 


62208 


+ 


-132840(3 iV°5? + 95093iV°5? + 206064(3iV25? + 863136(4A^^5? 
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+ 


+ 

+ 

+ 


1866240C5iV^5? - 20736387V25f - 2555712C3A^5f - m56UNgf 
233280C5Ngf - U28M0Ngf + 69984C3iV^c/[ff2 - 3 Q8AQN^glg2 
l9l28mC:iN‘^9lg2 - 69984C4A^^5l52 + 5598720C5A^^5l52 + 4710727V^5[52 
1664064C3iV5l52 + 1819584C4A^5l52 - f>fi^872QC^Nglg2 
38387UNglg2 - m2C:iN^g\gl + 3759N^gfgl + 445824(3N^gfg^ 
775555(4.N‘^g\gl + 29202587V2^?^^ _ i7infyQSC,3Nglgl 
225ll52CiNg\gl + 11664000C5A^5i5i - 6035150iVc/?5| 

454895(3N^g\gl + 209952C4iV'5f5i - m5254N^glgl - 754272(3Ng\gl 
4im94UNglgl - mmi59C^Ng\gl - l2927mNglgl + 14904C3A^'5^52 
85594N^g\gi - Q22533Q(3Ng\gi - 973944UNg\gi + 5531849C3Ng\gi 
14219955Ng\gl - 33 IUUC 3 Nglgl + Q29855 CaN glgl 
8398080C5iV5?52 + 893478iV5?ffi + 434159(3Nglgl - 189792C^Nglgl 
545952Nglgl - 4788072(35® + 332424(45^ + 1632960(55® 

l6^ + 


The immediate checks on these expressions are the internal ones similar to those alluded to 
for the wave function and coupling constant renormalizations. Specifically the higher order 
poles in e in the renormalization constants are not independent but depend on the lower order 
simple poles. That these higher order poles correctly emerge indicate that the procedure is not 
inconsistent. We note that unlike the expression given in (13.4p we have not included the wave 
function renormalization constants in the determination of the renormalization constants of the 
mixing matrix as is the usual procedure for a set of operators. See, for instance, m for a similar 
renormalization in QCD. 


A more appropriate check on our results rests in comparing with results from another quan¬ 
tum field theory. This is possible through the presence of the 0{N) symmetry which means that 
the anomalous dimensions can be extracted using another expansion method which is the large 
N expansion. While this was noted in m it is worth summarizing the background to the large 
N technique because of the overlap with the renormalization group at a critical point. If the 
/3-function has a non-trivial fixed point at the value 5 c, where this could also represent a vector 
of coupling constants such as we have here, then the renormalization group functions evaluated 
at 5 c are termed critical exponents which are renormalization group invariants. Moreover at 
a fixed point the critical exponents from field theories which are invariably different in nature 
can be the same. This universality is the key to our large N checks. Here the relevant scalar 
field theories which lie in the same universality class as the 0{N) nonlinear a model, 0{N) (j)‘^ 
theory and (|6.ip . The former two ordinarily reside in dimensions less than six and are respec¬ 
tively perturbatively renormalizable in two and four dimensions. In dimensions differing from 
their canonical dimension they may cease to be perturbatively renormalizable. Instead they are 
renormalizable above their canonical dimension in the sense that at their Wilson-Fisher fixed 
point in d-dimensions the critical exponents can be determined in the large N expansion. As N 
is a dimensionless parameter, which remains dimensionless in d-dimensions unlike the coupling 
constant in dimensionally regularized perturbation theory, then the quantity 1/N becomes a 
valid parameter for a perturbative expansion when N is large. Since the critical exponents can 
be computed to several orders in 1/A^ and in d-dimensions and because they are related to the 
critical renormalization group functions, they contain information on the anomalous dimensions 
of all theories in the same universality class at that fixed point. For the case we consider here, 
(Etd, the critical exponents corresponding to the two wave function renormalization constants 
and masses are known to 0(1/A^^) in d-dimensions and for <;//* to 0(1/A^®), [35l [36l [371 138]. 
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To be more specific and make contact with earlier work, [35l|36], it is worth recalling the situ¬ 
ation for 0{N) (j)'^ theory and define the relevant critical exponents. The form of the Lagrangian 
which is most appropriate is 

^ \ ( 6 . 11 ) 

where a is regarded as an auxiliary field which if eliminated produces the usual 4>^ Lagrangian. 
The single coupling constant g here appears with the quadratic term in a as it is in this particular 
form that the large N evaluation of the critical exponents is developed at high order in 1/N, 
[351 [Ml EZ]. Moreover, this formulation allows one to observe which theories he in the same 
universality class. For instance, the 0{N) nonlinear a model has a similar formulation but the 
final term is linear rather than quadratic in a and has a different coupling constant which is 
dimensionless in two dimensions. Moreover, there is a similarity to (|6.ip from the point of view 
of the interaction but the six dimensional theory has an additional coupling. However, only one 
of the critical points of (16.11) is in the same universality class as (16.lip . As a point of reference 
we use similar notation to [Ml [36] to define the full scaling dimensions of the fields. In [Ml |36| 
the critical exponents of and a were a and ^ respectively. It should be noted that these are 
not related to the physical exponents of earlier sections and are distinct to this section. We have 
modified the notation of [351IM] here for the full dimension of the fields to avoid confusion with 
the usual use of a and /3 as critical exponents. In terms of the anomalous contributions they 
are defined by 

a = — I + , P = 2 — 7] — X (6-12) 

where rj is the anomalous dimension of (/>* and y is the anomalous dimension of the interaction of 
(16.lip . The canonical dimension of j3 is in keeping with our mass operator canonical dimension 
convention. In terms of the renormalization group functions of (16.11) 

'y<p{gic, g 2 c) = , ia{gic,g 2 c) = - g - x ■ (6.13) 

Using ()6.4I) and (16.5p we have determined the location of fixed point in the large N expansion, 
{gic,g 2 c), and evaluated the field anomalous dimensions to the orders in 1/N to which they are 
known. It is satisfying to record that we find total agreement with our expressions. 

It is possible to repeat this check for Xij{gi,g 2 ) which requires some care. In the 4>‘^ formu¬ 
lation of the universal theory at the fixed point, (16.111) . the critical exponents of the two masses 
are straightforward to deduce. That for was discussed in early sections and in fact is equiv¬ 
alent to j3 above. The critical exponent for O 2 can be deduced from the final term of (I6.11|) . In 
relation to extracting information on the /3-function of 0(N} theory in the neighbourhood of 
four dimensions the relevant critical exponent is denoted by to and is proportional to the critical 
slope of the /3-function. However, in relation to the equivalent theory dehned in six dimensions 
that exponent would give the mass anomalous dimension of a as is evident from the form of the 
final term of (|6.1ip . Therefore, in terms of the exponents the critical point mass dimensions of 
the two mass operators are 


Ai = 2 - r] - X , A 2 = 2 w . (6.14) 

The remaining matter is to reconcile this argument with Jij{gi,g 2 ) at criticality and check if the 
exponents computed from it at the large N fixed point agree with the above known exponents 
at 0{1/N'^). The key to this is to compute the eigen-anomalous dimensions of 7ij(fl'i, 52 )- These 
are given by 

7 ±(fi'i, 5 ' 2 ) = 2 - [ 711 ( 51 , 52 )+ 722 ( 51 , 52 ) ±Arf( 5 ri, 52 )] (6.15) 
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where the discriminant is given by 


^d{gi,g 2 ) = [[711 ( 51 ,52) - 722 ( 51 ,52)]^ + 4712 ( 51 , ff 2)721 ( 51 ,52)] (6.16) 

and we have included the canonical dimension of 2 here. In [20] it appears that the canonical 
dimension of {d — 2) was used for the mass dimensions based on the dimensions of the con¬ 
stituent fields. However, here we retain the value of 2 as that convention is essential in getting 
consistency with this particular check between different theories. Evaluating 7 ±( 5 i,g 2 ) at the 
point {gic-,g 2 c) and expanding to 0{1/N‘^) we find that 7 -( 5 ic) 52 c) agrees exactly to four loops 
with Ai. Similarly 7_|_(g'ic, g' 2 c) is in precise agreement with A 2 to the same accuracy. In par¬ 
ticular the canonical dimension of 2 is crucial for ensuring the consistency of the latter as it 
derives from the exponent u in the 0{N) (j)'^ theory and that exponent near four dimensions 
corresponds to corrections to scaling rather than a mass operator. We regard these large N 
comparisons and in particular the second on the mass operators as non-trivial checks on our 
perturbative computations. It is worth mentioning that a one loop evaluation of a mass mixing 
matrix was given in [T8| as well as a two loop version for the related Sp{N) version of (16.ip in 
m- The renormalization group functions of the Sp{N) symmetric version of (16.Ih are related 
to those of the 0{N) theory by replacing N in each expression by {—N). However, as far as 
we can see both of the mass mixing matrix computations appear to have included 1-particle 
reducible diagrams. Low order in e checks of the critical exponents with respect to the large N 
mass operator exponents were discussed in [T8l [20] . 

One of the motivations of extending the three loop results of [20] to the next loop is to 
examine the critical point structure of the 0{N) two coupling theory. In [20] it was noted that 
there are various critical value of N for which the fixed point structure has different properties. 
For instance for N > Ncr there are three distinct critical points at real values of gi and <72 
which are perturbatively unitary. In addition for values of N in the range there are 

non-unitary fixed points. In |20] the values of these critical values of N were computed in an e 
expansion from the three loop /3-functions building on the one loop work of [251 [T8] . Therefore, 
we extend those estimates here using (|6.4p and (|6.5[) and use the same notation and method but 
in our conventions. First, we introduce the new scaled coupling variables x and y by 
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i 




(6.17) 


where we recall d = 6 — 2e. The presence of i here is to be consistent with [sniEi] given our 
coupling constant conventions. Then to find the critical couplings and Ncr one solves the set of 
equations 


/3i(5i,52) = /32(5i,52) = 


d/3i 8(32 d(3i 8132 


= 0 


(6.18) 


dgi 8 g 2 8 g 2 8 gi 

where the final equation is the condition for at least one zero eigenvalue of the Hessian of the 
matrix of derivatives of the /3-functions. Solving the resultant three equations perturbatively in 
d-dimensions there are three solutions which we designate A, B and C. The critical values in 
each of the three cases are 


= 1038.26605 - 1219.67959e - 1456.693326^ 3621.68482e^ O(e^) 

= 1.018036 - 0.01879e 0.0276066^ - 0.02587e^ O(e^) 

= 8.90305 - 0.42045e -k 4.067196^ - 2.00941e^ -k 0(6^^) (6.19) 


aJ = 1.02145 -k 0.06506e - 0.006526^ -k 0.20347e^ -k ©(e'^) 
x^c = *[0.23185 -k 0.17773e - 0.158226^ -k 0.61640e^ -k 0(6“^)] 

= /[0.25582 -k 0.22746e - 0.171066^ -k 0.77176e^ -k 0(6“^)] (6.20) 
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0(1) 

0(e) 

O(e^) 

0(6^) 

J. y 

1038.2660 

654.0820 

454.7593 

421.7574 

u^cr 

1.0180 

1.0087 

1.0126 

1.0133 

yA 

ycr 

8.9031 

8.6977 

8.8670 

9.5082 


Table 10. Estimates of critical value of N and location of fixed points for solution A in five 

dimensions using Fade approximants. 



0(1) 

0(e) 

O(e^) 

O(e^) 

cr 

1.0215 

1.0551 

1.0524 

1.0539 

^cr 

0.2318i 

— 

0.2933i 

0.3073i 

yB 

Ucr 

0.2558i 

- 

0.3385i 

0.3564i 


Table 11. Estimates of critical value of N and location of fixed points for solution B in five 

dimensions using Fade approximants. 


and 


AT^ = _ 0.08750 + 0.69453e - 3.53076e^ + 22.49021e^ + O(e^) 

= 0.13175 - 0.33427e + 0.48270e2 - 3.843496^ + O(e^) 

= - 0.03277 + 0.26911e - 1.43791e2 + 10.20700e^ + O(e^) . (6.21) 

We have checked that the results to O(e^) agree exactly with those given in [20]. With the 
O(e^) terms now present we can revisit the analysis of the location of the boundaries for the 
unitary versus non-unitary theories for five dimensions. In [20j a value of N^r = 64.253 emerged 
by setting e = 4 in N^. This is a significant distance from the leading value of 1038.266. As 
the O(e^) correction is both large and has a positive sign this means that the estimate for Wr 
will increase. Using the simple substitution approach then Wr = 516.963 which is larger than 
expectations from other methods. However, this large increase is perhaps more indicative of a 
lack of convergence and it seems more appropriate to estimate Wr by using Fade approximants. 
Using a [0, L — 1] Fade approximant for the Lth loop we find the values given in Table 10. 
It appears that the value of Wr settles to around 400 which while larger than the three loop 
estimate given in |20| is significantly smaller than the leading value. The corresponding values 
for solutions B and C are given respectively in Tables 11 and 12. In both instances the values 
of Ncr appear to converge to 1.05 and 0 respectively. Indeed for B the 0{e^) correction confirms 
the observation of [20] that > 1. 

One interesting application of our analysis is to compare estimates for the dimension of cj) 
with recent estimates using the conformal bootstrap method of [28]. There the N = 500 theory 
was considered directly in strictly five dimensions and the estimate of = 1.500409 was given. 



0(1) 

0(e) 

O(e^) 

O(e^) 

ly cr 

j^cr 

yC 

iJcr 

- 0.0875 
0.1317 

- 0.0328 

- 0.0176 

0.0581 

- 0.0064 

- 0.0082 

0.0347 

0.0039 

- 0.0035 

- 0.0112 

0.0228 


Table 12. Estimates of critical value of N and location of fixed points for solution C in five 

dimensions using Fade approximants. 
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This is remarkably close to the estimate obtained using the exponent evaluated to 0{1/N^) in 
the large N method which was = 1.500414. From the location of the /S-function zeroes for 
N = 500 we have 

X = 0.805458 + 0.276177e + 0.939812e2 + 8.0672426^ + O(e^) 
y = - 9.455850 - 4.724456e + 3.073550e2 + 140.713929e3 + O(e^) . (6.22) 

Using the four loop term and a Fade approximant we find the estimate of = 1.500537 

compared to a value of 1.500976 using the three loop expression. While not precisely on top of 
the other two methods there is a hint that a five loop computation may bridge the difference. 

As our final part of the analysis we briefly consider several low values of N in order to extend 
the results of [25l[Isl[2iTj where N = 0 and 1 were considered. The former value has been already 
considered in section 3 as it is the case of purely one field. For N = 1 there are two non-trivial 
fixed points one of which is stable at §2 = to leading order in e and the other is at gi = 52 
exactly, [251 [20] . For this unstable one the Lagrangian can be rewritten in such a way as to have 
a double copy of the basic theory (12.ip and it is believed that the flow is away from the unstable 
case to the infrared stable fixed point, pQj. Given our four loop analysis it is a simple exercise 
to extend the expressions for the critical exponents at each fixed point to the next order in e. 
At the g 2 = gi fixed point we have 

= 2 - l.lllllle - 0.1179706^ + 0.174760e^ - 0.631636e'^ + 0(e®) (6.23) 

for the dimensions of the fields and 

io+ = 2.000000e - 3.0864206^ + 12.725343e^ - 72.5220126"^ + 0(e®) 

Lo- = - 0.222222e - 0.2359406^ + 0.3495206^ - 1.263272e^ + 0(e®) (6.24) 

for the eigen-critical exponents of the matrix 

(6.25) 

9 i=gicg 2 = 92 c 

at the fixed point. In (|6.24ll we use the convention that d is not added to the eigenvalues of 

Uij and also note that the signs of the leading terms confirms the fixed point is unstable. In 

addition 

A+ = 2 - l.lllllle - 0.1179706^ + 0.174760e^ - 0.631636e^ + 0{e^) 

A- = 2 + 0.2222226 -b 0.235940e^ - 0.349520e^ + 1.263272e"^ -b 0(e®) (6.26) 

for the mass matrix eigen-exponents which shows A+ = A^-. By contrast at the infrared stable 
point we have 

A^ = 2 - 1.1002006 - 0.0937916^ -b 0.1605196^ - 0.545803e^ -b O(e^) 

A^ = 2 - 1.1222446 - 0.1435376^ -b 0.188846e^ - 0.721707e^ -b O(e^) 

UJ+ = 2.0000006 - 3.0927666^ -b 12.776556e^ - 72.867332e"^ -b O(e^) 
u- = 0.220441e -b 0.1750936^ - 0.316680e^ -b 1.2007816"^ -b 0(e®) 

A+ = 2 - 1.1222446 - 0.1435376^ -b 0.188846e^ - 0.7217076“^ -b 0(e®) 

A_ = 2 -b 0.100200e - 0.0464036^ - 0.156191e^ -b 0.336654e^ -b 0(e®) . (6.27) 

We note that our expressions for oji agree with m at three loops when the different conventions 
are accommodated. Clearly A_ = A^ as expected, |20|, in contrast to the situation at the 


UJid — 


dPi 

dgj 
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unstable fixed point. There the other eigen-exponent of the mass matrix was equivalent to the 
field dimension for = 1. In Appendix B we have recorded the field and mass eigen-exponents 
for a variety of values of N for the three fixed points that occur in the solution of f3i{gi,g2) = 0. 
The specific cases we considered are those chosen in [28] for a conformal bootstrap analysis. The 
same feature which has just been noted is apparent there. In other words at certain fixed points 
one of the mass eigen-exponents is equivalent to the anomalous dimension of the a field similar 
to the large N results. However, which mass exponent is identified depends on the nature of the 
underlying hxed point. 


7 Discussion. 

We close with various observations. First, we have carried out the four loop renormalization 
of theory in six dimensions. Clearly this has been a nontrivial exercise since we had to 
calculate the full set of master four loop massless 2-point functions to the requisite orders in e. 
Though tedious due to the large amount of integration by parts required, the master evaluation 
rested on the corresponding known masters in four dimensions, |l6|. That one can relate them 
was possible through the methods of mm and using general properties such as Weinberg’s 
theorem, m- One of the original aims was to refine the e expansion estimates of critical 
exponents for several physical problems. Overall the exponents we computed at the next order 
of precision are in reasonable accord with numerical approaches. Though the accuracy for low 
dimensions was driven by a deeper underlying property and that was the use of two dimensional 
conformal symmetry, m- There since critical exponents are known exactly one could constrain 
the exponent estimates and allow us to extract the behaviour across several dimensions. Such 
techniques now make the e-expansion reasonably competitive, but it is perhaps in the application 
to the more recent studies of 0{N) symmetric theories that will be beneficial in future. One 
issue examined in [20| was the range of N defining the conformal window which was suggested 
in pnj to drop from the one loop value of Wr = 1038, [25l |T8|, to around Ncr = 64. Here 
we took a more conservative approach in applying summation techniques to suggest that while 
the bounding value drops, it maybe does not reduce so far as this. While successive three and 
four loop estimates suggest a value settling to around Ncr = 400 it would be premature to 
regard this problem as having been resolved. It may be that the perturbative approach is not 
as fully equipped to give a definitive answer in comparison to, say, the conformal bootstrap 
machinery. Also as noted in m the N = 1 theory may be related to the deformed (3,10) 
minimal conformal field theory in two dimensions. Therefore, if true then one could perform a 
similar analysis across the gap to six dimensions. It would be interesting then to try and match 
predictions with physical systems. 

One intriguing possibility is the potential application of the universality between 2 and 6 
dimensions to non-scalar theories such as those with gauge symmetry or supersymmetry. Indeed 
the conformal window of QCD is of interest, [Hj. Akin to the connections of four dimensional cf)'^ 
theory and six dimensional theory one question which would be worth considering in fnture 
is what if any is the higher dimensional theory which is in the same universality class as the 
Banks-Zaks fixed point in QCD. To an extent there is already a parallel in QCD with what has 
been discussed here and in m for the large N connection in the scalar field theories. At the 
Wilson-Fisher fixed point in 2 < d < 4 it is known that QCD is in the same universality class 
as the non-abelian Thirring model in the large Nf expansion, [68]. Here Nf is the number of 
(massless) quarks and it is important to appreciate that the connection is with respect to this 
particular parameter being large as opposed to Nc being large where Nc is the number of colours 
in SU{Nc) QCD. Like (j)'^ theory QCD is perturbatively renormalizable in four dimensions and 
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the non-abelian Thirring model is renormalizable in two dimensions. Indeed this fact has been 
used to compute large Nj critical exponents to varying orders in l/iVf in the Thirring model and 
demonstrate that their e expansion with respect to four dimensions agrees with the analogous 
renormalization group functions of QCD. See, for example, iMiiin]. While parallel to the lower 
end of the chain of nonlinear a model, (j)'^ and (j)^ theories which are renormalizable in two, four 
and six dimensions respectively, what is missing in the QCD instance is the corresponding six 
dimensional theory. To extend the QCD chain would require an understanding of the critical 
point properties of a spin-1 field in contrast to a spin-0 field. We have referred to this field not 
as a gauge held or gluon because in the two dimensional non-abelian Thirring model the spin-1 
held plays the role of an auxiliary held rather than a gauge held much in the same way that 
one replaces the interaction in (j)^ theory by a 3-point vertex at the expense of introducing an 
auxiliary held. In (16.1|] this appears as a and is the connecting held for the chain in 4 < d < 6 
as well as being the Lagrange multipler held in the two dimensional 0{N) nonlinear a model. 
There it imposes the constraint that the helds lie on a multidimensional sphere. So the a held 
is effectively the lynchpin held which also underlies the chain of theories at the Wilson-Fisher 
hxed point across the dimensions. 
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A Master Integrals. 


In this appendix we record the explicit values of the master integrals required for the four loop 
renormalization to various orders in e near six dimensions. Terms beyond the 0(l/e) are needed 
due to the presence of spurious poles in e from the integration by parts algorithm of Laporta, 
[42] . We have not included all the basic integrals used but only those which are not products of 
lower loop integrals nor which contain only simple self-energy subgraphs in order to save space. 
These are straightforward to construct directly by expanding products and ratios of Euler F- 
functions. We use the same notation used in [46] for the definition of the topology. The explicit 
graphs of the integrals are given in [46|. We have 
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(A.l) 


In several cases we have given terms up to the level of (^7 despite the fact that this number does 
not appear in the final renormalization group functions. It would be present in the finite parts 
of the various Green’s functions. For compactness we have expressed the master integrals in 
(lA.lIl in the corresponding G-scheme in six dimensions which is the standard way to present 
them, [46]. In other words for each loop order the factor 


G = 


r(l + e)(r(l-e))- 
r(2-2e) 


(A.2) 


is included. This circumvents the need to include, for instance, the Euler-Mascheroni constant 7 
which would otherwise appear throughout and which is ordinarily included in part in the change 
from the minimal subtraction to modified minimal subtraction scheme, m- 


B Field dimensions at various fixed points. 


In this appendix we give the various critical exponents for (16.11) to O(e^) at the same fixed points 
value of N as given in the appendix of [28]. This is partly for completeness but also to note 
interesting structure. For the various values of N considered there are three distinct fixed points 
which are labelled in the same way as [28] . These were termed ‘critical’, ‘theory 2’ and ‘theory 3’ 
and we retain that nomenclature here but label the exponents with the respective superscripts 
c, 1 and 2. The first case considered in [28] was N = 600 and the three sets of exponents are 

= 2 + (0.000036i - 0.996357)e + (0.000142i - 0.008332)e2 

+ (-0.000852i + 0.001801)e^ + (-0.002328i + 0.014101)e'^ + 0(e®) 
a;; = 2 + (-0.060308i + 0.167456)e + (0.172529i- 0.263035)e2 

+ (0.274055i + 0.000784)e^ + (-0.846566i + 0.464339)e'^ + 0{e^) 

A% = 2 +(-1.594272 - 0.673433i)e +(-0.231855 + 1.280485i)e2 

+ (0.890531 + 3.353446i)e^ + (1.853548 - 1.733472^)6"^ + O(e^) 

Ai = 2 +(0.167456 - 0.060308i)e +(-0.263035 + 0.172529^)6^ 

+ (0.000784 + 0.274055i)6^ + (0.464339 - 0.846567i)6"^ + 0{e^) (B.l) 

A^ = 2 + (-0.000036i - 0.996357)6 + (-0.000142i - 0.008332)6^ 

+ (0.000853i + 0.001801)6^ + (0.002328i + 0.014101)6"^ + 0(6®) 

Al = 2 + (0.060308i + 0.167456)6 + (-0.172529^- 0.263035)6^ 

+ (-0.274055i + 0.000784)6^ + (0.846566i + 0.464339)6^ + 0(6®) 

A^ = 2 + (-1.594272 + 0.673433^)6 +(-0.231855 - 1.280485i)62 

+ (0.890531 - 3.353446^)6® + (1.853548 + 1.733472i)6"^ + 0(6®) 

A?. = 2 +(0.167457+ 0.060308i)6 +(-0.263035 - 0.172530^)6^ 

+ (0.000784 - 0.274055^)6® + (0.464339 + 0.846566i)6"^ + 0(6®) (B.2) 
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and 


Al = 2 - 0.997795e - 0.00219e2 + 0.003555e^ + 0.034922e^ + 0(e®) 

= 2 - 0.159643e + 0.167857e2 + 0.763443e^ + 4.639369e^ + 0(e®) 

Aij. = 2 - 0.159642e + 0.167857e2 + 0.763443e^ + 4.639369e^ + O(e^) 

Ai = 2 + 0.632464e + 0.506339e2 + 4.963931e^ + 35.781014e'^ + 0(e®) . (B.3) 

Here and for the expressions we give for the other values of N we note that 

= 7±(5'1c,5'2c) (B.4) 

includes the canonical dimension and these mass exponents definitions are the same as those 
used for the large N comparison with the large N critical exponents. For this value of N and 
others two of the points have critical exponents which are complex conjugates. Next the parallel 
values for N = 1000 are 

A^ = 2 + (-O.OOOOlOi - 0.997921)e + (0.000167i - 0.004378)e2 

+ (0.001518i + 0.000406)6^ + (0.021951z + 0.007247)e^ + O(e^) 

A^ = 2+ (0.010791z +0.115748)6 + (-0.186373z-0.172907)e2 

+ (-1.392417Z - 0.116672)6^ + (-21.819202i - 0.017897)6"^ + 0{e^) 

A^ = 2 +(-1.424029 + 0.158680z)6 +(-0.198690 - 2.582892^)6^ 

+ (-0.142920 - 22.664631z) 6^ + (-2.551351 - 349.159818z)6"^ + 0(6®) 

Ai = 2 +(0.115748 + 0.01079H)6 +(-0.172907 -0.186373 z)62 

+ (-0.116672 - 1.392417i)6^ + (-0.017897 - 21.819201i)6"^ + 0(6®) (B.5) 

Al = 2 + (O.OOOOlOi - 0.997921)6 + (-0.000167z - 0.0043778)6^ 

+ (-0.001518Z + 0.000406)6^ + (-0.021951z + 0.007247)6^ + 0(6®) 

Al = 2+ (-0.01079H +0.115748)6 + (0.186373^-0.172907)6^ 

+ (1.392417i - 0.116671)6® + (21.81920H - 0.017897)6^ + 0(6®) 

A^ = 2 +(-1.424029 -0.158680z)6 +(-0.198690 + 2.582892^)6^ 

+ (-0.142920 + 22.664631i)6® +(-2.551351 + 349.159818z)6^ + 0(6®) 

A?_ = 2 +(0.115748 - 0.01079H)6 +(-0.172907+ 0.186373i)62 

+ (-0.116671 + 1.392417i)6® + (-0.017897 + 21.819202i)6"^ + 0(6®) (B.6) 

and 

= 2 - 0.9986126 - 0.0015356^ + 0.0018646® + 0.0189096^ + 0(6®) 

A® = 2 - 0.1278486 + 0.1382896^ + 0.5984526® + 3.4999636^ + 0(6®) 

A®_ = 2 - 0.1278486 + 0.1382896^ + 0.5984526® + 3.4999636^ + 0(6®) 

Ai = 2 + 0.5256206 + 0.3862596^ + 3.8060346® + 27.0889096"^ + 0(6®) . (B.7) 

Finally, for N = 1400 we find 

Al = 2 - 0.9985266 - 0.0029246^ - 0.0000566® + 0.0040146^ + 0(6®) 

A^ = 2 + 0.0677156 - 0.1431316^ + 0.0053746® + 0.3484396^ + 0(6®) 

A)^ = 2 - 1.7650856 - 0.7674836^ + 0.2914376® + 2.7281776^ + 0(6®) 

Al = 2 + 0.0677156 - 0.1431316^ + 0.0053746® + 0.3484396^ + 0(6®) (B.8) 
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Al = 2 - 0.998575e - 0.002939e2 + 0.000342e^ + 0.0058776"^ + 0{e^) 

Al = 2 + 0.115801e - 0.123622e2 - 0.2790432e3 - 0.776356e^ + O(e^) 

= 2 -0.924175e + 0.427891e2- 1.113372e3-10.101218e'^ + 0(e®) 

A‘i = 2 + 0.115801e - 0.1236220e2 - 0.279043e^ - 0.776356e^ + 0(e®) (B.9) 

and 

= 2 - 0.998982e-0.001191e2 + 0.001224e3+0.012696e^ + 

A^ = 2 - 0.109951e + 0.120821e2+ 0.507233e^+ 2.895587e'^ + 

= 2 - 0.109951e + 0.120821e^ + 0.507233e^ + 2.895587e'^ + 0(e®) 

Ai = 2 + 0.462254e + 0.321631e^+ 3.184495e^+ 22.495823e^ + O(e^) . (B.IO) 

We note that the expression for A^ corrects an obvious typographical error in equation (A.23) 
of [28]. There the 0(e) term is not recorded although its actual coefficient appears as the 
coefficient of the O(e^) term. For the N = 1400 set the three fixed points again produces 
real critical exponents. Unfortunately in each set improving the series convergence using Fade 
approximants only applies to A^ along the lines discussed in the main text for N = 500. In 
each case we have recorded the mass mixing matrix eigen-critical exponents in our conventions 
as there is an interesting feature which extends the observation in the large N comparison in 
section 6. In each set of exponents and values of N A^- is equivalent to one of the mass eigen- 
critical exponents. In other words for finite values of N the field critical exponent and its mass 
exponent are equivalent. However, the particular eigen-exponent the field dimension equates 
to depends on the specific fixed point. In each of the cases presented here A^- corresponds to 
the minus exponent for the points designated critical and theory 2 but to the plus exponent for 
theory 3. Of the three only theory 3 has real exponents and this picture tallies with the large 
N checks discussed earlier. 
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